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ABSTRACT 


This  is  a  two-part  technical  report  that  considers  the  short-time 
statistical  modeling  of  atmospheric  noise  at  VLF/LF  frequencies  and 
applies  a  new  noise  model  to  the  performance  evaluation  of  various 
generic  VLF/LF  communication  modems.  In  Part  One  atmospheric  noise  is 
examined  mainly  from  a  statistical  point  of  view.  Next,  various  statis¬ 
tical  models  for  this  type  of  noise  are  examined;  finally,  a  new  model 
is  developed  and  compared  with  available  dat:..  In  Part  Two  this  noise 
model  is  applied  to  the  evaluation  of  the  error  probability  calculations 
of  various  generic  VLF/LF  modems.  Much  of  Part  Two,  particularly  in 
regard  to  a  receiver  system  employing  no  clipping,  is  based  on  classic  1 
white  Gaussian  noise  analysis  applied  to  generic  modems  of  interest. 

When  nonlinear  clipping  is  applied,  the  detailed  statistical  behavior  of 
atmospheric  noise  becomes  more  important,  and  results  depend  on  the  de¬ 
tailed  structure  of  the  noise. 


til 


CONTENTS 


ABSTRACT .  iii 

LIST  OF  ILLUSTRATIONS .  vii 

PART  ONE 

ATMOSPHERIC  NOISE  MODEL 

I  INTRODUCTION  .  3 

II  PHYSICAL  DESCRIPTION  OF  ATMOSPHERIC  NOISE .  5 

A.  Introduction .  5 

B.  The  Signals  from  a  Single  Lightning  Flash  .  6 

1 .  The  Flash  to  Earth .  6 

2.  The  Intracloud  Discharge  .  9 

C.  Characteristics  of  a  Thunderstorm  .  9 

D.  Climatology  of  Thunderstorms .  10 

E .  Summary .  10 

III  SUMMARY  OF  EXISTING  MODE! A5  FOR  ATMOSPHERIC  NOISE  ....  13 

IV  LOG-NORMAL  NOISE  MODEL  .  17 

A,  Introduction .  17 

B,  Properties  of  the  Log-Normal  Model .  20 

1.  First-Order  Statistics . 23 

2.  Higher-Order  Statistics .  31 

C,  Conclusion .  36 

PART  TWO 

ERROR  PROBABILITIES  OF  GENERIC  VLF/LF  MODEMS 

I  INTRODUCTION  .  39 

II  PROPERTIES  OF  THE  NOISE  MODEL .  41 

III  SIGNAL  REPRESENTATION .  47 

IV  COHERENT  BINAR\  RECEIVER  MODEL  .  49 


v 


V  PROBABILITY -OF  -ERROR  CALCULATIONS . 

A.  Binary  On-Off  Keying  (CW) .  55 

L.  Binary  Frequency-Shift  Keying  (FSK)  .  61 

C. -  Binary  Phase-Shift  Keying  (PSK) .  70 

D.  Binary  Minimum-Shift  Keying  ulSK)  .  76 

E.  M-ary  Minimum  Shift  Keying . 86 

VI  CONCLUSION .  93 

APPENDIX .  95 

REFERENCES .  99 

DISTRIBUTION  LIST .  103 

DD  Form  1473 


vi 


ILLUSTRATIONS 


Figure  1  Field  Changes  at  100  km  Due  to  Typical 

Lightning  Discharges .  7 

Figure  2  Comparison  of  Spectra  at  100  km  Due  to 

Various  Processes  in  a  Lightning  Flash.  ...  8 

Figure  3  Diurnal  Variation  of  Lightning  Occurrences.  .  11 

Figure  4  Envelope  Model  for  Gaussian  and  Log-Normal 

Processes .  19 

Figure  5  Probability  Distribution  of  the  Envelope  for 

Average  Voltage  2:  800  nV/m .  25 

Figure  6  Probability  Distribution  of  the  Envelope  for 

Average  Voltage  2:  300  and  50  pV/m .  26 

Figure  7  Envelope  Level  Crossings  and  Associated  Step 

and  Counting  Function  Processes  .  27 

Figure  8  Average  Number  of  Times  Envelope  Crosses 

Level  B .  32 

Figure  9  Probability  Distribution  of  Time  Between 

Crossings  of  Level  B.  .  .  .  . .  35 

Figure  10  Coherent  Receiver  Model — No  Clipping .  49 

Figure  11  Coherent  Receiver  Model-Clipping .  49 

Figure  12  Noncoherent  Detection  of  CW  Pulse — No 

Clipping .  56 

Figure  13  Binary  Error  Probability  for  On-Off  Xeying — 

No  Cl  ipping .  59 

Figure  14  Noncoherent  Detection  of  CW  Puise--Cl ipping  .  60 

Figure  15  Binary  Error  Probability  for  On-Off  Keying — 

Clipping .  62 

Figure  16  Coherent  Detection  of  FSK — No  Clipping.  ...  63 

Figure  17  Binary  Error  Probability  for  Coherent  FSK — 

Clipping .  66 


vii 


Figure  18  Noncoherent  Detection  of  FSK — No  Clipping  .  .  67 

Figure  19  Noncoherent  FSK  Detector  Representation  ...  67 

Figure  20  Binary  Error  Probability  for  Noncoherent 

FSK— Clipping .  71 

Figure  21  =  Binary  Error  Probability  for  Coherent 

PSK— No  Clipping .  73 

Figure  22  MSK  Modulator  Representation . 76 

Figure  23  Coherent  MSK  Detector — No  Clipping.  ......  77 

Figure  24  Binary  Error  Probability  for  Coherent  MSK — 

No  Clipping .  78 

Figure  25  Binary  Error  Probability  for  MSK--Clipping .  .  79 

Figure  26  Binary  Error  Probability  for  Differential 

MSK .  86 

Character  Error  Probability  of  MSK  Sequences.  89 


Figure  27 


PART  ONE 

ATMOSPHERIC  NOISE  MODEL 


1 


PART  ONE 


ATMOSPHERIC  NOISE  MODEL 

I  INTRODUCTION 

In  Part  One  we  develop  an  analytical  model  for  atmospheric  noise 
that  is  characterized  by  its  ’'impulsive''  phenomena.  Specifically,  a 
new  model  is  developed  for  radio  noise  that  originates  in  lightning 
discharges  and  propagates  large  distances  to  VLF/LF  receiving  systems. 
The  emphasis  here  is  on  development  of  an  accurate  short-time,  statis¬ 
tical  noise  model  that  can  be  easily  applied  to  the  evaluation  of  the 
performance  of  existing  VLF/LF  communication  systems,  so  that  measured 
noise  parameters  (such  as  average  power  and  deviation  of  the  noise 
voltage  envelope)  and  signal  parameters  can  adequately  specify  system 

performance.  Noise  parameters  and  estimation  of  noise  parameters  are 

,  * 

discussed  in  the  final  report  on  this  project.1  Our  goal  here  is  to 
model  the  short-time  statistical  characteristics  of  atmospheric  noise  as 
accurately  as  is  necessary  for  evaluating  the  performance  of  VLF/LF 
communication  systems  of  interest. 

The  modeling  approach  used  here  is  based  on  the  observation  that  a 
Gaussian  noise  model  is  inappropriate  for  received  atmospheric  noise, 
primarily  because  Gaussian  noise  does  not  have  the  large  dynamic  range 
exhibited  by  the  received  noise  envelope.  Generally,  the  received 
atmospheric  noise  envelope  has  a  log-normal  distribution  for  large 
envelope  values  and  a  Rayleigh  distribution  for  smaller  dynamic  ranges. 
Since  the  envelope  of  a  Gaussian  process  is  Rayleigh  distributed,  atmo¬ 
spheric  noise  has  a  Gaussian  behavior  only  for  small  dynamic  ranges. 

In  a  communication  system  limited  by  atmospheric  noise,  the  large  noise 
envelopes  have  the  most  influence  on  system  performance,  and  so  it  is 
important  to  characterize  accurately  the  large  dynamic  range  behavior  of 
atmospheric  noise. 

• 

References  are  listed  at  the  end  of  this  report. 
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We  next  discuss  the  background  for  the  noise  modeling  problem.  In 
particular,  we  examine  the  physical  characteristics  of  atmospheric  noise. 
Then  we  discuss  existing  models  for  atmospheric  noise  and  introduce  the 
noise  model  developed  in  this  study.  Finally,  the  properties  of  this 
new  noise  model  are  examined  and  compared  with  available  atmospheric 
noise  data. 


II  PHYSICAL  DESCRIPTION  OF  ATMOSPHERIC  NOISE* 

A .  Introduction 

The  main  source  of  atmospheric  noise  at  frequencies  in  the  VLF 
(3-30  kHz)  and  LF  (30-300  kHz)  bands  is  the  lightning  discharge.  Each 
flash  generates  a  s.  iuence  of  VLF  pulses  having  certain  characteristics. 
The  thunderstorm — as  a  whole--has  a  lifetime  during  which  the  nature  of 
the  lightning  flashes  that  take  place  and  their  rate  of  production  may 
vary.  The  occurrence  of  thunderstorms  has  a  climatologically  predictable 
spatial  and  temporal  pattern.  Thus  the  definition  of  the  thunderstorm 
sources  of  VLF  atmospheric  noise  may  be  considered  under  three  divisions: 
the  signals  due  to  a  single  flash,  the  sequence  of  discharges  in  a  single 
thunderstorm,  and  the  distribution  temporally  and  spatially  of  the 
centers  of  thunderstorm  activity  over  the  globe. 

After  the  VLF  noise  signals  have  been  generated,  they  are  propagated 
in  accordance  with  the  established  laws  of  VLF  radio  propagation.  Thus 
in  order  to  determine  the  noise  environment  at  a  particular  time  and 
locality,  three  steps  are  necessary:  first,  to  establish  the  nature 
of  the  noise  signals  emanating  from  a  typical  thunderstorm  center;  next, 
to  determine  the  relative  activities  and  the  distances  of  the  main 
thunderstorm  centers  from  the  locality  being  considered  at  the  specified 
time;  and,  finally,  to  introduce  the  modifications  caused  by  propagation 
to  the  noise  originating  in  the  various  thunderstorm  centers. 

This  chapter  summarizes  some  of  the  aspects  of  VLF  atmospheric 
noise  specified  above.  Particular  attention  is  given  to  source  effects. 
Most  of  the  topics  summarized  here  are  considered  in  much  greater  detail 
in  Ref  .  2  . 


This  chapter  was  prepared  by  Dr.  E.  T.  Pierce. 
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B. 


The  Signals  from  a  Single  Lightning  Flash 


1 .  The  Flash  to  Earth 

The  duration  of  a  discharge  to  earth  varies  from  a  few  tens 
of  milliseconds  to  over  a  second;  400  ms  is  perhaps  a  reasonable  mean 
value.  The  flash  is  characterized  by  an  initial  L  stage  during  which 
a  leader  moves  downward  from  the  cloud  towards  the  earth.  This  is 
succeeded,  after  contact  has  been  made  with  the  ground,  by  an  inter¬ 
mediate  stage  characterized  by  the  occurrence  of  one  or  more  brilliant 
return  (R)  strokes.  After  the  last  return  stroke  there  is  a  further 
gradual  flow  of  current  during  the  final  (F)  stage.  The  leader  stage 
of  the  flash  to  earth  occupies  some  50  ms,  while  100  ms  is  a  typical 
duration  for  the  F  stage;  thus  st  of  the  discharge  is  occupied  by  the 
intermediate,  R,  stage. 

Figure  1  illustrates  the  structure  of  the  VLF  noise  due  to  a 
flash  to  earth  at  a  distance  of  100  km.  The  field-change  records  in  the 
figure  [(a)]  are  typical  of  those  obtained  using  a  bandwidth  of  1-1000 
Hz;  these  records  are  dominated  by  the  electrostatic  component  of  the 
change  in  field.  The  records  of  Type  (b)  are  for  a  bandwidth  of  1-100 
kHz.  It  is  easily  seen,  from  comparing  records  of  Type  (a)  and  Type 
(b),  that  the  main  VLF  pulses  are  associated  with  return  strokes  (R 
pulses),  but  that  small  L  pulses  occur  during  the  leader  stage  while 
somewhat  larger  K  pulses  are  present  in  the  intermediate  and  final 
stages . 

The  number  of  return  strokes  per  flash  to  earth  may  vary  from 
one  to  ten  or  more;  the  average  is  about  three  to  four.  The  median 
interval  between  the  strokes  is  some  50  ms,  and  the  distribution  of 
intervals  is  skew  towards  the  large-interval  end;  some  15  percent  of 
the  intervals  exceed  100  ms,  while  about  20  percent  are  less  than  30  ms. 
Figure  2  illustrates  the  typical  amplitude  spectrum  for  an  R  pulse;  the 
vertical  scale  applies  at  a  distance  of  100  km.  Note  that  the  VLF  dis¬ 
turbance  due  to  the  first  return  stroke  is  larger  in  magnitude  and  peaks 
at  a  higher  frequency  than  that  caused  by  a  subsequent  return  stroke. 
However,  the  sizes  of  both  the  first  and  the  subsequent  R  pulses  vary 
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FIGURE  2  COMPARISON  OF  SPECTRA  AT  100  km  DUE  TO  VARIOUS 
PROCESSES  IN  A  LIGHTNING  FLASH 


over  a  wide  nnge;  if  the  amplitudes  are  expressed  in  decibels  relative 
to  the  median,  the  distribution  is  normal  with  a  standard  deviation  of 
about  6  dB . 

The  K  pulses  occur  intermittently  throughout  the  intermediate 
and  final  stages.  The  average  time  interval  between  the  pulses  is 
about  10-15  ms,  and  the  distribution  of  time  intervals  is  skew  towards 
the  large  value  end.  A  typical  flash  to  earth  produces  perhaps  25  K 
pulses.  The  spectrum  of  an  average  K  pulse  tends  to  peak  at  a  rather 
higher  frequency — see  Fig.  2 — than  that  of  an  R  pulse,  while  the 
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magnitude  is  perhaps  some  10  percent  of  that  for  an  R  pulse.  The  ampli¬ 
tude  distributions  for  the  K  and  R  pulses  appear  to  be  similar,  both 
being  log-normal  with  a  standard  deviation  of  some  6  dB.  Thus,  since 
there  are  perhaps  ten  times  as  many  K  pulses  as  R  pulses  per  flash,  it 
is  not  unusual  for  the  largest  K  pulse  in  a  discharge  to  ground  to  be 
larger  than  the  smallest  R  pulse. 

Irregular  L  pulses  occur  during  the  L  stage  of  the  discharge. 
Typically  two  or  three  pulse  trains,  each  consisting  of  some  50  L  pulses 
with  a  total  duration  of  about  3  ms,  are  produced,  '.’he  disturbance 
associated  with  an  individual  L  pulse  is  small,  and  since  the  L  pulses 
are  generated  for  a  time  normally  less  than  10  ms,  the  L  effects  can 
usually  be  neglected  by  comparison  with  the  K  and  R  phenomena. 

2  .  The  Intracloud  Discharge 

The  structure  of  the  VLF  noise  from  a  flash  within  a  thunder¬ 
cloud  is  much  simpler  than  that  of  the  disturbance  associated  with  a 
discharge  to  earth.  Intracloud  flashes  do  not  generate  R  pulses,  but 
they  produce  intermittent  K  pulses  throughout  the  duration  of  the  dis¬ 
charge.  In  addition,  small  irregular  pulses  are  generated  in  the  initial 
stages  of  the  intracloud  flash,  but  their  effect  is  negligible  compared 
with  those  of  the  K  pulses. 

The  durations  of  flashes  to  earth  and  those  of  intracloud  dis¬ 
charges  are  comparable.  Furthermore,  the  time  separation  and  magnitudes 
of  the  K  pulses  ere  similar  for  the  two  types  of  flash.  Thus,  to  a  first 
approximation,  we  may  regard  the  VLF  noise  due  to  an  intracloud  discharge 
as  being  identical  to  the  contribution  for  a  flash  to  earth  associated 
with  K  pulses . 

C .  Characteristics  of  a  Thunderstorm 

An  average  thunderstorm  has  a  duration  of  perhaps  2  hours.  During 
the  lifetime  of  the  storm,  the  flashing  rate  varies  appreciably;  the 
peak  rate  is  about  ten  discharges  per  minute  and  the  average  some  three 
per  minute.  Most  of  the  discharges  in  a  storm  are  of  the  intracloud  or 
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allied  types;  typically,  some  20  percent  of  all  flashes  go  to  earth,  and 
this  proportion  increases  with  an  increase  in  geographical  latitude. 

D .  Climatology  of  Thunderstorms 

It  is  usually  estimated  that  about  2000  thunderstorms  are  active 
at  any  instant  and  that  the  rate  of  lightning  occurrence  integrated  over 
the  whole  surface  of  the  earth  is  about  100  per  second.  These  figures 
are  exactly  compatible  with  an  average  flashing  rate  per  storm  of  three 
per  minute . 

There  are  three  especially  important  global  centers  of  thunderstorm 
activity:  in  Southeast  Asia,  Central  Africa,  and  South  America.  Each 

center  has  its  maximum  of  activity  in  the  local  afternoon;  this — see 
Fig.  3 — leads  to  a  characteristic  global  diurnal  variation  in  lightning 
occurrence.  This  variation  has  of  course,  superimposed  seasonal  changes. 

E .  Summary 

The  procedures  for  estimating  the  atmospheric  noise  at  a  particular 
locality  and  time,  due  to  global  lightning,  are  now  clear.  We  first 
determine  the  rate  of  flashing  at  each  of  the  main  thunderstorm  centers 
from  curves  such  as  those  of  Fig.  3.  If  for  any  center  the  flashing 
rate  is  5r  per  second,  then  we  have  about  4r  intracloud  flashes  and  r 
discharges  to  earth.  Each  flash  to  earth  gives — on  an  average — three 
to  four  R  pulses,  separated  by  some  50  ms,  while  all  discharges,  whether 
to  earth  or  not,  generate  about  25  K  pulses  at  intervals  of  perhaps  10 
ms.  We  may  assume  the  flash  occurrence  within  a  main  thunderstorm 
center  to  be  random,  since  many  individual  storms  will  be  active  simul¬ 
taneously.  Thus  at  peak  activity  one  of  the  main  thunderstorm  centers 
will  be  producing  rbcut  5\-  flashes  per  second,  and  therefore  some  35  R 
pulses  and  perhaps  1250  K  pulses  per  second.  The  bigger  R  pulses  will 
give  the  log-normal  high-amplitude  end  often  apparent  in  an  experi¬ 
mentally  observed  amplitude  distribution  curve;  the  smaller  R  pulses 
and  the  K  pulses  will  yield  the  smal 1 -ampl itude  Rayleigh  end  of  the 
amplitude-distribution  curve. 
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FIGURE  3  OIURNAL  VARIATION  OF  LIGHTNING  OCCURRENCES 

After  the  characteristics  (amplitude  and  temporal  variations  in 
pulses  emitted)  of  a  particular  source  have  been  established,  the  contri 
bution  of  this  source  to  the  noise  at  the  locality  being  considered  can 
be  derived  by  introducing  the  appropriate  propagational  modification  to 
the  source  disturbance.  Integration  over  all  sources  then  gives  the 
total  noise.  Note  that,  if  all  the  significant  source  centers  are  very 
distant  or  strongly  active  at  intermediate  distances,  then  the  received 
noise  will  be  dominated  by  a  very  large  incidence  of  smallish  pulses. 
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Only  if  the  activity  is  principally  due  to  local  storms,  or  to  medium 
activity  at  intermediate  distances  will  the  incidence  of  large  isolated 
pulses  be  significant. 
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Ill  SUMMARY  OF  EXISTING  MODELS  FOR  ATMOSPHERIC  NOISE 


Before  a  new  model  for  atmospheric  noise  is  described,  a  brief 
summary  of  available  models  is  now  presented.  Various  analytical  models 
for  received  atmospheric  noise  have  been  proposed,  usually  from  one  of 
two  points  of  view.  Briefly,  the  most  interesting  approach  from  a 
physical  point  of  view  takes  the  received  noise  to  be  the  weighted  sum 
of  contributions  from  individual  lightning  discharges.  Although  this 
approach  is  physically  well  justifiable,  it  has  the  disadvantage  that 
the  resulting  models  are  not  analytically  tractable.  The  alternative 
approach  is  an  empirical  method  that  yields  analytical  models  chosen  to 
fit  the  measured  statistics  of  the  noise.  Unfortunately,  not  only  is 
direct  physical  support  for  this  procedure  lacking,  but  the  existing 
empirical  models  suffer  from  the  fact  that  only  the  first-order  statis¬ 
tics  of  the  noise  are  considered  while  the  higher-order  statistics  are 
neglected . 

.The  models  based  on  the  weighted  sum  of  individual  lightning  dis¬ 
charges  are  often  referred  to  as  f iltered-impulse  models. 3-6  These 
physically  motivated  noise  models  have  the  typical  form 

a(t)  3 i  P\X  '  ti)  >  (1) 

i 

where  the  ja.^  ate  independent,  identically  distributed  random  variables 
whose  distribution  is  deduced  from  lightning  discharge  statistics  and 
propagation  consideration;  p(t)  gives  the  form  of  the  noise  pulse  re¬ 
sulting  from  an  individual  lightning  stroke,  as  shaped  by  the  front-end 
filter  of  the  receiver;  and  the  jt^J  ore  the  occurrence  times  of  the 
individual  lightning  strokes.  Under  the  assumption  that  each  individual 
lightning  stroke  is  independent  of  others,  various  results4-*  have  been 
obtained  on  the  first-order  statistics  of  the  noise.  In  particular, 
probability  distributions  arc  typically  derived  for  those  models,  since 
this  distribution  has  been  measured  quite  extensively.4*0-10  Although 
the  analytical  results  due  to  different  assumpCions  on  p(t),  |t  |,  and 
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|  a^|  result  in  varying  agreement  with  measured  data,  the  resulting 
analytical  probability  distributions  are  quite  complicated  and  difficult 
to  use  in  performance  calculations.  Indeed,  the  resulting  probability 
distributions  often  cannot  be  written  in  closed  form  except  for  some 
limiting  cases.  Also  the  assumption  of  independence  between  \ncTlvidual 
lightning  strokes  which  is  crucial  to  the  solution  of  the  first-order 
statistics  of  a(t)  does  not  appear  to  be  true  at  VLF 7  This  means  that 
any  second-order  (or  higher-order)  statistic  would  not  be  accurately 
modeled  by  these  noise  models.  Conceivably,  the  independent  assumption 
might  be  eliminated  by  introducing  dependence  in  the  f iltered-impulse 
model,  as  was  done  by  Furutsu  and  Ishida,4  who  considered  the  case  of 
Poisson-Poisson  noise.  Although  this  interdependence  model  has  not  yet 
been  demonstrated  to  yield  the  distribution  of  interlevel-crossing 
intervals  observed  by  Watt  and  Maxwell  at  VLF,7  it  is  already  too  diffi¬ 
cult  to  handle  for  performance  calculations. 

Noise  models  based  on  the  empirical  method  generally  consist  of 
mathematical  expressions  constructed  to  fit  the  measured  ta  on  the 
first-order  statistics  of  the  envelope  of  the  received  noise.  These  are 
models  of  the  envelope  distribution  and  not  of  the  whole  noise  process 
as  is  the  case  with  f iltered-impulse  models.  Not  only  are  these  em¬ 
pirical  models  not  phy  ^ically  motivated,  but  by  their  construction  they 
are  good  for  only  first-order  calculations.  Although  they  are  relatively 
simple  to  handle  for  first-order  calculations,  their  range  of  application 
is  limited. 

Ibukun*~  presents  a  good  summary  of  various  empirical  models  that 
have  been  proposed.  Several  other  workers* * ~l 3  have  had  some  success 
in  finding  physical  justification  for  empirical  models.  Perhaps  the 
most  important  of  these  models,  because  of  its  simplicity  and  its  close¬ 
ness  of  fit.  to  measured  data,  is  the  model  that  takes  the  envelope  of 
tho  received  atmospheric  noise  to  be  Rayleigh  distributed  at  low  values 
of  the  envelope  and  log-normally  distributed  at  high  values.  Beckmann* 3 
has  given  a  good  physical  argument  which  supports  this  model,  particu¬ 
larly  in  the  situation  where  there  is  little  local  thunderstorm  activity. 
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Several  workers10 >'  3 j1 4  have  proposed  models  similar  to  that  considered 
by  Beckmann,  although  they  differ  in  regard  to  how  much  of  the  two  dis¬ 
tributions  should  be  combined  to  result  in  the  best  model. 

One  of  the  most  important  recent  atmospheric  noise  models  is  due  to 
Hall15  and  is  called  the  generalized  "t"  model.  This  modeling  effort 
takes  a  point  of  view  different  than  those  leading  to  either  the  filtered- 
impulse  models  or  the  empirical  models  discussed  above.  The  generalized 
"t"  model  describes  the  received  atmospheric  noise  as 

a(t)  .  b(t)n(t)  ,  (2) 

where  n(t)  is  a  narrow-band  Gaussian  process,  and  b(t)  is  a  slowly 

varying  process,  independent  of  n(t).  This  model  gives  some  physical 

support  to  perhaps  the  simplest  of  the  empirical  models,10'16  which  has 

probability  P  (V)  that  the  noise  envelope  exceeds  the  value  V, 
o 

Vv)  -  [‘  Mr)']  -  <3) 

where  V  is  the  average  value  of  the  envelope  and  3  and  r  are  two 
parameters  to  be  chosen.  This  empirical  model  is  also  the  same  as  the 
model  studied  by  Mertz.*">lJ  Not  only  do  the  first-order  statistics  of 
Hall's  generalized  "t"  model  agree  closely  with  this  class  of  empirical 
models,  but  it  has  the  advantage  over  empirical  models  in  that  it  can 
be  specified  to  give  a  good  fit  to  the  higher-order  statistics  of  the 
noise  also. 

I'sing  the  approach  of  Hall,  wc  next  develop  an  atmospheric  noise 
model  that  is  perhaps  mathematically  simpler  than  Hall's  generalized 
“t"  model  and  yet  has  as  good  an  agreement  with  first-  and  second-order 
statistical  data  as  the  generalized  ”t .” 


15 


IV  LOG-NORMAL  NOISE  MODEL 


A .  Introduction 

Turning  attention  now  to  the  problem  of  developing  a  model  for  re¬ 
ceived  atmospheric  noise,  we  note  that  this  noise  is  always  observed 
through  the  passband  of  some  receiver  filter.  Now,  if  the  receiver  is 
sufficiently  narrow-band,  the  noise  at  the  receiver  output  can  reasonably 
be  assumed  to  be  modeled  well  as  a  Gaussian  process.  This  follows  from 
the  fact  that  narrow-band  filtered  noise  is  the  sum  of  contributions  from 
many  independent  lightning  discharges,  none  of  which  is  dominant  at  the 
filter  output.  Experimental  data  indicate,  however,  that  the  bandwidth 
required  to  achieve  this  condition  at  VLF  is  less  than  50  Hz,  so  a 
Gaussian  assumption  is  not  always  physically  viable  at  VHF.  The  goal  of 
the  model  development  here  is  the  formulation  of  an  analytical  model 
that  is  both  an  accurate  description  of  the  received  noise  and  suitable 
for  application  to  the  calculation  of  VLF/LF  communication  system  per¬ 
formance.  As  far  as  the  communication  problem  is  concerned,  it  appears 
necessary  to  model  the  atmospheric  noise  prior  to  any  receiver  opera¬ 
tions,  so  a  Gaussian  assumption  is  not  always  justified.  This  is  par¬ 
ticularly  true  when  nonlinear  operations  are  performed  by  the  receivers. 

The  modeling  problem  can  be  simplified  by  noting  that  for  communica¬ 
tion  application  the  receiver  bandwidths  are  substantially  smaller  than 
the  band  center  frequency.  This  fact  enables  the  received  atmospheric 
noise  to  be  regarded  as  a  narrow-band  random  process.  This  is  always 
satisfied  for  communication  problems  and  is  not  nearly  as  strong  an 
assumption  as  a  Gaussian  assumption.  Almost  all  the  available  experi¬ 
mental  data4  have  been  obtained  in  narrow-band  conditions. 

Measured  data  on  atmospheric  noise  indicate  thi t  this  noise  has  a 
Gaussian  behavior  at  low  amplitudes  and  a  log-normally  distributed 
envelope  for  large  amplitudes.  Beckmann  explains  this  fact  by  noting 
that  measured  atmospheric  noise  usually  consists  of  the  effects  of  many 
lightning  discharges  around  the  world.  When  no  single  discharge  domi¬ 
nates  at  any  instant  of  time,  then,  apj lying  the  central  limit  theorem, 
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we  should  expect  a  Gaussian  behavior.  On  the  other  hand,  when  a  par¬ 
ticular  individual  discharge  dominates,  the  measured  amplitude  has  the 
statistical  characteristics  of  the  individual  discharge  which  is  essen¬ 
tially  log-normal  in  character.  Since  the  larger  amplitudes  have  the 
greater  influence  on  the  performance  of  any  communication  system,  our 
model  emphasizes  the  log-normal  characteristic  of  atmospheric  noise. 

In  particular,  we  model  atmospheric  noise  as  a  narrow-band  process  with 
a  log-normal  envelope  with  the  form 

a(t)  =  A  e”^  sin|u)ot  +  Q ( t ) J  ,  (4) 

where  n(t)  is  a  stationary  Gaussian  process  with  zero  mean  end  auto¬ 
correlation  given  by 


R  (T)  -  n(t)n(t  +  r)  ,*  (5) 

n  ’ 

A  is  a  constant  (to  be  determined  from  noise  power  estimates),  and  Q(t) 
is  a  random  phase  process  independent  of  the  Gaussian  process  n(t) . 

To  illustrate  the  behavior  of  the  envelope  of  this  model,  consider 

Fig.  4,  In  Fig.  4(a),  we  plot  a  sample  Gaussian  process  n(t)  .  In  Figs. 

n  ( t ) 

4(b),  (c),  and  (d)  we  plot  e  to  show  how  the  envelope  of  this  model 

behaves  for  typical  values  of  (voltage  deviation,  the  ratio  of  the 

root -mean-square  voltage  to  the  average  voltage  of  the  noise  envelope) . 
By  proper  adjustment  of  parameters,  this  model  reflects  the  impulsive 
nature  of  atmospheric  noise  that  has  passed  through  a  receiver  front- 
end  filter. 

In  our  model  given  by  Eq .  (4)  we  have  not  yet  specified  the  be¬ 
havior  of  the  phase  process  q(t).  Measurements  of  the  instantaneous 
frequency  distribution  for  atmospheric  noise  show  that  it  is  similar  to 
the  frequency  distribution  of  narrow-band  Gaussian  noise.  Hence,  oven 


*  _ 

Here  x  is  the  expected  value  of  x. 
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TIME  -  orbitrory  unit* 

0’0«)«5 


FIGURE  4  ENVELOPE  MODEL  FOR  GAUSSIAN  AND  LOG  -  NORMAL 
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thoi  ,h  the  envelope  behavior  of  atmospheric  noise  and  that  of  Gaussian 
noise  differ  considerably  for  large  envelope  values,  the  phase  and  fre¬ 
quency  behavior  are  quite  similar.  We  therefore  assume  that  © ( t )  be¬ 
haves  like  the  phase  of  narrow-band  Gaussian  processes.  This  means  that 
at  any  given  instant  t  ,  the  phase  0(to)  is  a  random  variable  that  is 
uniformly  distributed  over  [0,  2rr  1  and  is  independent  of  the  envelope. 

B .  Properties  of  the  Log-Normal  Model 

The  log-normal  narrow-band  noise  model  for  atmospheric  noise  is 
given  by 

a(t)  =  A  en^  sinlx  t  + 

L  0 

Assuming  stationarity  we  note  that  the  envelope  of  this  noise, 

E(t)  =  A  en(t)  ,  (7) 

has  autocorrelation  given  by 

R^t)  =  E(t)E(t  +  y) 

=  A2  en(t)+n(t+T)  (8) 

„  r o2+r  (t)*i 

*  A2  eL  °  n  J  , 

2 

where  a  -  R  (0)  is  the  variance  of  n(t)  . 
n  n 

The  average  power  of  the  noise  model  is  given  by 
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R  (0) 
a 


2  2n(t) 

=  A  e  s 


2  2°n  1 

=  A  e  X  g 


2  r 

in  uu  t  +  Q  ( t ) 
° 


,2  2q 
A  n 
“e 


(10) 


If  we  assume  that  the  noise  has  an  approximately  flat  energy  distribution 

of  N  watts  per  hertz  across  the  receiver  front-end  bandwidth  of  W  Hz, 
o 

then  we  have  the  relation 


2  2a2 

N  W  =  e  " 
o  2 


(ID 


Consider  the  average  and  rms  voltages  of  the  envelope.  The  average 


is 


E  =  E ( t)  =  A  e 
ave 


n  (t) 


=  A  e 


B  n 


(12) 


and  the  rms  envelope  voltage  is 


rms 


(E2(t) 


A  e 


|a2  J2n(t)£ 

) 


-  '  A 

-  )A 
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(13) 


The  voltage  deviation,  V. ,  is  defined  as 

d 


|Erms  | 


Vd  ■  20  lo»10lE— I 
i  ave  i 


=  20  lo*io’e 


l.i’nl 

I  I 


(14) 


=  10rin  1Og10  6 
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2 

Knowledge  of  V.  gives  us  the  value  of  a  .  Similarly,  knowing  V,  (hence 
2  d  n  d 

o  )  and  the  noise  power  density,  N^,  we  can  find  the  parameter  A.  Values 

of  V,  and  of  F  ,  which  is  directly  related  to  N  ,  are  estimated  in  CCIR 
d  a'  o’ 

reports  ,19 

We  have  now  related  parameters  of  our  noise  model  directly  to  the 
measurement  parameters  that  are  typically  used  in  world  CCIR  maps.  Our 
next  task  is  to  show  how  well  this  noise  model  agrees  with  available 
statistical  data  on  atmospheric  noise.  The  measurements  required  to 
check  the  validity  of  the  model  fall  into  two  categories. 

The  first  of  these  categories,  and  the  one  for  which  the  greatest 
amount  of  experimental  data  is  available,  is  concerned  with  the  first- 
order  statistics  of  the  random  process.  The  particular  measurements 
in  this  category  that  have  been  reported  in  the  literature  are  measure¬ 
ments  of  the  probability  distribution  of  the  envelope  of  the  received 
noise,4?7-10  and  measurements  of  the  average  number  of  level  crossings 
per  unit  time  of  a  fixed  level  by  the  envelope  of  the  received  noise.7?9 
This  latter  calculation  has  not  been  carried  out  for  any  of  the  empirical 
models  discussed  earlier,  although  Nakai20  has  obtained  numerical  re¬ 
sults  in  agreement  with  the  experimental  data  of  Watt  and  Maxwell  for  a 
f iltered-impulse  model  in  which  the  noise  pulses  occur  in  a  Poisson 
manner.  Hall15  has  also  obtained  good  agreement  with  experimental  data 
with  his  generalized  "t"  noise  model.  In  fact,  he  has  also  shown  good 
agreement  with  second-order  experimental  data  which  we  discuss  next. 

The  second  category  of  measured  data,  and  the  one  where  far  fewer 
data  are  available,  is  concerned  with  the  second-  and  higher-order 
statistics  of  the  random  process.  The  particular  measurements  that  have 
been  reported  in  this  category  are  measurements  of  the  probability  dis¬ 
tribution  of  the  time  interval  between  crossings  of  a  specified  level  by 
the  envelope  of  the  noise.7  While  available  experimental  data  are  sparse, 
being  restricted  to  a  few  measurements  at  VLF,  it  is  also  true  that  the 
analytical  derivation  of  these  statistics  is  complicated,  requiring 
machine  computation  in  the  general  case.  Hall  was  able  to  compare  a 
limiting  case  of  his  generalized  "t"  noise  model  with  second-order  data 
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and  demonstrated  reasonable  agreement  between  his  model  and  the  available 
data.  We  use  similar  approx  imat  ioi^s  ,  to  demonstrate  agreement  of  our 
log-normal  model  with  higher-order  data. 

1 .  First-Order  Statistics 

The  envelope  of  our  log-normal  narrow-band  atmospheric  noise 
model  is  given  by 


E(t) 


=  A 


n(t) 

e 


(15) 


We  now  derive  some  first-order  statistical  properties  of  this  envelope 
and  compare  them  with  available  data.  At  any  instant  of  time, 

E  =  A  e"  (16) 


is  a  random  variable,  where  n  is  a  zero-mean,  Gaussian,  random  variable 
2 

with  variance  afl.  The  probability  distribution  of  E  is 


P(V)  =  P  fE  £  Vi 
r 

>  PrjA  en  5  V| 

=  Pr{en  <;  (V/A)} 


=  Pr{n  ^  loge(V/A)| 


=  4 


(log  (V/A) 
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T" 


where 


P  [A]  is  the  probability  of  the  event  A,  and 
r  2 

-a  _•!. 

2 

ijf(a>  =  I  ^5=  e  dx 
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is  the  well-known  distribution  of  a  zero-mean,  unit-variance  Gaussian 
random  variable.31  Defining 


P  (V)  =  1  -  P(V) 
o 


=  1 


flog@(V/A) 


V? 


(19) 


we  compare  this  probability  distribution  of  the  envelope  derived  from 
our  model  with  measured  data  on  the  distribution  of  the  envelope  of  at¬ 
mospheric  noise.  Figures  5  and  6  compare  the  above  distribution  with 
experimental  data  for  various  receiver  parameters.  Note  that 
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V  =  E  =  A  e 
ave  ave 


(20) 


Hence 


Vv>  =  *1 


riog  (v/v  |  -  Jo 

e\  ave'  2  n 

T 


(21) 


and  a  is  the  only  free  parameter  in  the  figures, 
n 

The  plotted  results  are  reasonably  self-explanatory,  showing 
good  agreement  between  the  model  and  the  measured  data.  Again,  we  point 
out  that  our  intention  is  to  model  the  large  envelope  variations  more 
accurately  than  previously,  since  these  variations  have  the  most  in¬ 
fluence  on  communication  system  performance.  Consequently,  our  model 
matches  the  experimental  data  more  accurately  for  larger  values  of  the 
envelope . 

We  next  consider  another  experimentally  measured  first-order 
statistical  property,  the  average  frequency  of  envelope  level  crossings. 
Consider  the  envelope  process  and  the  level  B  indicated  in  Fig.  7(a). 

We  are  concerned  with  the  average  number  of  times  E(t)  crosses  the  level 
B  in  T  seconds.  This  is  equivalent  to  asking  for  the  average  number  of 
times  the  Gaussian  process  n(t)  crosses  the  level 


FIGURE  5  PROBABILITY  DISTRIBUTION  OF  THE  ENVELOPE  FOR  AVERAGE  VOLTAGE  s  800  pVI m 


FIGURE  7  ENVELOPE  LEVEL  CROSSINGS  AND  ASSOCIATED  STEP 
AND  COUNTING  FUNCTION  PROCESSES 


b  =  log  (B/A) 
e 


In  Fig,  7(b)  the  equivalent  Gaussian  process  crossing  level  b  is  shown. 


Wc  define  a  counting  function 


c  ( t )  =  |  u  [n ( t)  -  bj ! 


I  u  [£(  t)  -  B]l 


which  gives  a  unit  impulse  each  time  [n(t)  ^  b  or  E(t)  =  B]  the  level  b 
is  crossed  by  n(t) .  Since  we  want  to  count  the  number  of  crossings  in 
a  time  interval  of  length  T,  the  count  is  given  by 


H  (T)  =  /  c(t)  dt 


|u[n(t)  —  b ] I  dt 


Using  the  well-known  relationship 


we  get 


VT) 


I  6[n 

Jo 


(t)  -  b]  I n ( t ) I  dt 


Since  we  have  a  stationary  Gaussian  random  \  nriabl  n(t)  and  n(t)  are 
statistically  independent.  This  follows  from  the  relation 


n  ( t )  h  ( t )  - 


1  d 

2  dT  n  (t) 


)  d  2  . 

2  dT  n  (t) 


0 


Hence, 


N  (T)  =  6  L n ( t )  -  b]  ri( t) 1  T 
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(28) 


The  random  variable  n(t)  has  a  density  function 
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The  probability  density  of  ri(t)  is  found  by  noting  that  the  derivative 
is  a  linear  operation  and  that  linear  operations  on  a  Gaussian  process 
result  in  another  Gaussian  pr*_~ess.  Thus  n(t)  is  a  Gaussian  random 
variable  with  zero  mean  and  variance  determined  as  follows: 


R(-)  -  n(t)n(t  +  r) 


(31) 


R(-)  --  V  R(r) 
d  t 


(32) 


-  n  ( t )  n  (  t  +  r ) 


Bv  stationaritv  we  have 
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so 


R(r)  =  R(t) 
dr 


(34) 


n(t  -  r)n(t) 


Hence 


R( 0)  =  -  [n( t)  ] 

and  the  probability  density  of  ii(t)  is 
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(35) 


p.U)  = 


_  \_  /  1 

s„,o>]  i  2[-s„l0>]i 


(36) 


The  term  Iri(t)  |  is  giver,  by 


|ri(t)  |=2/  -v  P^('v)  do 
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(37) 


do 


Since  R  ( 0)  is  related  to  how  rapidly  the  envelope  process  goes  through 
n  _ 

changes,  the  quantity  jri(t)j  depends  on  higher-order  statistics  of  the 
envelope  process . 

Combining  Eqs.  (28)  and  (30)  we  find  that  the  average  number 
of  times  the  envelope  of  our  noise  model  crosses  the  level  B  in  T 
seconds  is  given  by 
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For  the  special  case  of  T  -  1,  wo  have 

wm  |  [ 

*■" *  ■*'  exp  -  ^ 

2-r-  ( 

n  ' 


nb(1) 


log  (B/A) 
e 


2  - 


111 


(39) 


30 


Figure  8  compares  this  average  level  crossing  with  experimental  data. 

It  shows  a  basic  agreement  of  the  results  derived  by  using  our  log-normal 
model  with  measured  level-crossing  results.  Again,  we  note  that  the 
agreement  is  somewhat  better  for  the  larger  values  of  the  envelope  in¬ 
tensity.  The  significance  of  this  agreement  is,  of  course,  that  it  is 
further  verification  of  the  applicability  of  the  model  ac  far  as  the 
first-order  statistics  of  the  noise  envelope  are  concerned.  From  a 
physical  point  of  view,  this  agreement  increases  confidence  in  the  log¬ 
normal  narrow-band  noise  model,  since  it  indicates  that  on  the  average 
the  envelope  of  the  received  noise  is  accurately  modeled.  To  complete 
the  verification  of  the  applicability  of  this  noise  model  for  atmospheric 
noise,  it  remains  to  investigate  the  manner  in  which  these  variations 
occur  with  time.  This  is,  of  course,  a  function  of  the  second-  nnd 
higher-order  statistics  of  the  noise,  wnich  is  the  next  topic  of 
d iscussion . 

2 .  Higher-Order  Statistics 

Completion  of  the  verification  of  the  applicability  of  the 
narrow-band  log-normal  process  to  received  atmospheric  noise  requires 
investigation  into  the  higher-order  statistics  of  the  noise.  Although 
the  average  rate  of  level  crossings  of  the  envelope  of  our(  model 
generally  agrees  with  observed  atmospheric  noise,  we  must  yet  verify 
that  the  higher-order  statistics  of  the  model  can  be  so  chosen  that  the 
relationship  between  the  noise  process  at  various  distinct  times,  as 
predicted  by  the  model,  is  consistent  with  measured  results.  As  men¬ 
tioned  earlier,  the  available  experimental  data’  that  are  dependent  on 
the  higher-order  statistics  of  the  noise  consist  of  measurements  of  the 
probability  distribution  function  of  the  interval  between  successive 
crossings  of  a  fixed  level  by  the  envelope  of  the  noise.  Inspection  of 
these  data  indicates  that  at  VLF  the  noise  pulses  do  not  in  general 
occur  in  a  Poisson  fashion  but  rather  that  received  noise  pulses  are 
usually  statistically  dependent  on  preceding  ones.  Generally,  calcu¬ 
lation  of  the  probability  distribution  of  the  interlevel-crossing  time 
interval  is  difficult,  requiring  numerical  techniques  to  obtain  even  an 
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FIGURE  8  AVERAGE  NUMBER  OF  TIMES  ENVELOPE  CROSSES  LEVEL  B 


approximate  solution . 22-24  This  is  true  even  for  the  simplest  case  of 
Gaussian  processes,  which  is  the  only  case  that  has  been  treated  in  any 
detail  in  the  literature ,a3>24  Even  in  the  Gaussian  case,  the  problem 
is  further  restricted  by  consideration  of  Markov  processes  and  a  few 
other  specific  spectra.  Therefore,  we  will  not  attempt  to  find  an  exact 
solution  for  the  process  of  interest  here  but  will  resort  to  simplifying 
assumptions  in  order  to  obtain  an  approximate  solution. 

Consider  the  envelope,  E(t),  crossing  some  level  B,  and  let  T 
be  defined  as  the  interval  between  a  down-crossing  at  the  level  B  and 
the  next  up-crossing  of  the  same  level.  (This  is  the  statistical  data 
measured  by  Watt  and  Maxwell.7)  The  probability  that  T  exceeds  Tq  is 
given  by 


pit  i  =  p  It  >  t  [ 

o'  o  '  r »  o’ 

=  Pp|E(t)  <;  B,te[0,To]|E(0-)  >  B,E(0+)  b}  (40) 

Pr{E(t)  <;  B,te[0,ToJ;E(0-)  >  B,E(0+)  <  b} 

=  Pr|E(0-)  >  B,E(0+)  <  B  |  5 


where  E(0+)  is  the  value  of  E(t)  at  a  small  increment  of  time  after 
t  *  0,  and  E(O-)  is  the  corresponding  value  at  a  small  increment  of 
t ime  before  t  =  0  . 


At  this  point,  we  make  several  simplifying  assumptions  by 
considering  N  time  samples  of  the  envelope  process  and  using  these  time 
samples  to  represent  the  process.  Taking 


V  V  tN 


(41) 


as  the  N  +  1  time  samples,  we  approximate  the  desired  probability  by 
PjEltjl  *  B,E(  t2  )  S  B,...,E(tN)  S  B,E(to)  >  B} 


P  j  T  ) 
o  o ' 


Pr|Eil0>  >  B;EI  Y!  *  B| 


.  ( 42) 
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Recall  that  our  envelope  process 


E(t)  =  A  e 


n(t) 


(43) 


is  band-limited  to  W,  so  in  T  seconds  there  are  approximately 

’  o 


N  =  2W  T 


(44) 


independent  samples.  Taking  samples  to  be  independent,  we  then  have 

N 


PrlE(‘ol  ^  E(  7T  Pr)E|V  5  *4 


P  (T  ] 
o'  o 1 


k=l 


P.J]E|tJ  >  B 


k=2 


From  Eq .  (17)  we  have 


E|tJ  <i  B 


J]  prlE(V  s  B! 


(45) 


P  |Ei 
r* 


B}  =  ’]r 


(log  (B/A) 


therefore 


P  ( T  )  2: 
o'  o' 


(2WT  -1) 
o 


(46) 


(47) 


In  Fig.  9  we  plot  this  distribution  along  with  data  taken  from 
Watt  and  Maxwell.  The  assumption  of  independent  samples  neglects  certain 
dependencies  between  adjacent  level  crossings.  This  assumption  was 
necessary  to  obtain  a  simple  closed-form  expression  for  Pq(Tq)  ancl 
not  necessarily  a  limitation  of  the  log-normal  model.  Indeed,  by  appro¬ 
priately  choosing  R^t),  a  more  accurate  agreement  with  data  is  possible, 
and  only  our  ability  to  compute  Pq(To)  in  a  convenient  form  is  limited. 

In  a  more  general  case,  for  example,  we  would  have 
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FIGURE  9  PROBABILITY  DISTRIBUTION  Or  TIME  BETWEEN  CROSSINGS  OF  LEVEL  B 


* 


p 

o  o 


(48) 


where 

b  -  loge(B/A) 

and 


K 


l  it. 
n  \  1 


is  an  (N  +  1) -d imensional  covariance  matrix. 


C .  Conclusion 

We  conclude  that  the  log-normal  atmospheric  noise  model  gives  good 
agreement  with  measured  data,  especially  for  larger  values  of  the  enve¬ 
lope.  This  applies  for  the  first-order  statistics  as  well  as  the  higher- 
order  ones.  This  model  is  also  relatively  easy  to  handle  analytically 
and  thus  is  useful  for  calculating  the  performance  of  communication 
systems . 
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PART  TWO 


ERROR  PROBABILITIES  OF  GENERIC  VLF/LF  MODEMS 


PART  TWO 


ERROR  PROBABILITIES  OF  GENERIC  VLF/LF  MODEMS 

I  INTRODUCTION 

In  this  part  we  discuss  the  performance  of  various  generic  VLF/LF 
communication  systems  of  interest.  The  usual  performance  measure  of 
probability  of  error  will  be  expressed  in  terms  of  the  signal  strength 
and  noise  parameters  at  the  receiver's  antenna.  As  we  have  discussed 
in  Part  One,  the  noise  has  an  impulsive  non-Gaussian  behavior  which  we 
have  modeled.  The  receiving  systems  we  will  consider  can  also  have 
nonlinear  clipping  which,  together  with  the  non-Gaussian  noise  behavior, 
makes  computation  of  the  error  probabilities  difficult  if  not  impossible 
without  some  assumptions.  Our  assumptions  are  discussed  and  justified 
as  they  arise  in  the  development  of  error  probability  expressions.  We 
first  discuss  the  noise  model,  the  signal  characteristics,  and  the  re¬ 
ceiving  system  models  before  carrying  out  the  development  of  the  error 
probability  expressions.  The  VLF/LF  communication  systems  we  consider 
include  binary  frequency  shift  keying  (FSK),  phase  shift  keying  (PSK), 
and  minimum  shift  keying  (MSK) ,  all  operating  in  a  coherent  mode.  On- 
off  keying  (CW)  and  noncoherent  FSK  are  also  considered,  as  well  as 
differential  phase  shift  keying  (DPSK)  .  Two  signaling  schemes  that  use 
long  sequences  of  MSK  pulses  (chips)  are  also  discussed. 


39 


II  PROPERTIES  OF  THE  NOISE  MODEL 


At  VLF  and  LF,  the  dominant  noise  source  is  generally  atmospherics. 
Local  activity  often  produces  more  noise,  but  this  is  much  more  variable 
and  unpredictable.  In  our  application,  receivers  will  usually  be  located 
in  areas  far  from  man-made  noise  sources  and  high  thunderstorm  areas. 

We  therefore  assume  that  the  noise  is  due  to  atmospherics,  which  are  due 
to  distant  thunderstorm  activity. 

The  short-time  noise  model  we  use  assumes  the  form 


a(t) 


=  A  e 


n(  t) 


cos^nf^t  +  9(t)j 


(49) 


where  n(t)  is  a  stationary  Gaussian  process,  9 ( t )  is  a  stationary  random 
process  independent  of  n(t),  and  A  is  a  constant.  As  we  ha«e  shown  in 
Part  One,  this  noise  model  agrees  well  with  available  data,  including 
second-order  level-crossing  data.  The  agreement  is  particularly  good 
for  large  fluctuations  of  amplitude.  Since  large  fluctuations  of  noise 
amplitude  are  more  important  than  average  noise  in  terms  of  causing 
errors  in  any  system,  it  is  particularly  important  to  model  this  region 
accurately.  We  note  that  for  small  fluctuations  the  experimental  data 
show  the  noise  to  behave  like  a  Gaussian  process.  Although  our  model 
does  not  give  Gaussian  r.oise  at  small  fluctuations,  wt  '.nut  use  of  the 
above  experimental  result  as  a  justification  for  or.e  of  rur  assumptions. 

Consider  next  some  of  the  properties  of  this  noise  model  and  their 
relation  to  measured  noise  parameters.  The  phase  random  variable  at  time 
t  is  uniformly  distributed  over  {0,  2").  Hence  at  any  time  we  have 


!  2  2nd)  2T  T 

(t)  -  A  c  '  cos  2-fQt  ♦  ?(t> 


2  1  n 
A  —  c 


(50) 


■n 


c 


where 


cos^j"  Srrf^t  +  0(t) 

2  T~7 

a  =  n  (t) 


(51) 

(52) 


Although  we  have  assumed  stationr.rity  in  our  short-time  noise  model  ,  it 

is  clear  that  there  are  diurnal  and  seasonal  variations  in  the  noise 

characteristics.  These  slow  variations  can  be  accounted  for  in  the  two 

2 

basic  parameters  of  our  short-time  noise  model  given  by  A  and  . 

n 

2 

To  relate  A  and  o  to  important  measured  noise  parameters,  we  now 
note  that  the  noise  model  is  valid  for  the  atmospheric  noise  that  has 
been  filtered  through  a  receiver  front-end  filter  of  bandwidth  W.  In 
our  application,  W  is  between  200  Hz  and  1000  Hz.  Assuming  that  the 
atmospheric  noise  oefore  filtering  is  approximately  flat  in  spectrum 
over  the  filter  bandwidth,  we  have 

_  2,2 

a2(t)  =  ~  A2  o  n  =  jN'  W  (53) 

2  o 


where  is  the  noise  spectral  density  per  unit  bandwidth.  The  envelope 
of  the  process  is 


n  ( O 

E(t)  =  A  e  (54) 

which  has  an  average  value 

i.,2 

E(t)  ,  A  e  n  (55) 


and  an  rms  value 


-  a 


2 

n 

n 

e 


( 56) 


A  measure  of  the  variability  of  the  noise  fluctuations  is  the  quantity 

V,,  defined  as 
d7 
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If  the  envelope  E(t)  is  restricted  to  be  less  than  or  equal  to  B.  then 
the  resultant  envelope  c(t)  is 


c(t) 


A  e^^  if  not  H(B) 


if  H(B) 


(63) 


with  probability  density 


(x)  =  5(x  -  B)Pb  +  j^l  -  ug(x)j  p^log^x/A)  j 


(64) 


where 


exp!-  ^~2  (  > 
y 2"-n  I  2"n  I 


p  (x)  -  - exP 


Ur(x)  = 


1  x  >  B 


0  x  <  B 


(65) 


(66) 


and  5(x)  is  the  Dirac  delta  function.  The  mean  square  of  this  "clipped" 
envelope  is 

log  (B/A) 


L 


>2  (,  /  |  2 
■jp  exp  3a  -  — — ,  da  +  B  P 


.2  n 
=  A  e 


A  e 


2  log  (B/A) 
2a  I  e 


7^1  r^l 

.  |(-2=n)2) 

l  7^rp  (^n 

1  /log  (B/A)  -  2  -2  \  ) 


da  +  R  PQ  (67) 
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Using  Eq .  (61)  for  P  we  have 

B 


_  9  2 

2  2 
c  =  A  e 


1  -  $ 


floge(B/A)  -  2on, 


2 

+  B  $1 


floge(B/A)’ 


(68) 


This  noise  model  is  used  to  compute  the  desired  error  probabilities, 
2 

The  parameters  A  and  q  are  obtained  irom  measured  and  predicted  values 

of  F  and  V,.  Equation  (68)  is  used  to  approximate  the  effects  of 
ad 

clipping  on  the  received  noise. 
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Ill  SIGNAL  REPRESENTATION 


The  modulation  system  to  be  used  at  VLF  is  strongly  limited  by  the 
characteristics  of  the  high-power  amplifiers  and  high -efficiency  antennas 
that  constitute  present  stationary  VLF  transmitters.  The  most  desirable 
modulation  signals  would  maximize  the  information  transmitted  in  the 
allowed  bandwidth  through  a  noisy  channel  for  given  fixed  peak  power. 

Given  a  peak  power  limitation,  the  highest  average  power  is  radiated 
when  the  transmitted  wave  has  a  constant  envelope.  An  additional  factor 
that  necessitates  a  constant  envelope  is  the  final  RF  amplifier.  Existing 
high-power  VLF  facilities  use  Class  C  amplifiers  because  of  their  effi¬ 
ciency;  this  mode  of  operation  requires  constant  envelope  and  continuous 
phase  for  highest  efficiency.  Discontinuity  in  envelope  or  phase  of  an 
RF  wave  also  causes  large  transients  in  the  transmitter  output.  Tran¬ 
sient  response  of  the  transmitting  system  depends  upon  the  system  band¬ 
width,  which  at  VLF  is  generally  limited  by  the  high  Q  inherent  to  the 
most  efficient  transmitting  antennas. 

Because  of  practical  limitations  on  the  transmitters  at  VLF  and  LF, 
we  consider  only  angle-modulated  signals,  with  emphasis  on  those  that 
are  ph  ase  continuous.  In  particular,  we  restrict  our  attention  to 
binary  FSK,  PSK,  and  MSK  signaling  schemes.  One  exception  that  we  con¬ 
sider  is  on-off  CW  keying. 

During  a  signaling  time  of  T  seconds,  one  of  two  signals  will  be 
transmitted.  Lotting  Hq  be  the  hypothesis  that  a  signal  corresponding 
to  "zero"  is  transmitted  and  bo  the  hypothesis  that  a  "one"  signal  is 

sent,  we  have 


as 


H  : 

o 


the  two  types  of 


m  (t)  r.  ”2p  cost  2rf  t 
o  L  o 

m.  (t)  -  COS  2rrf  t 

i  L  O 

signal.  For  FSK  we  have 


*  e0(»] 

*  ¥'>] 


(89) 

(70) 
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and 


(71) 


e  (t)  =  Au)t 
o 


Q1  <  t)  =  -  /\uot 

Note  that  the  average  power  for  both  signals  is  P,  and  let  p  be  the 
correlation  between  m_(t)  and  m1(t).  Hence 

m  ( t)m, (t)  dt  .  (72) 

o  1 

The  signals  we  consider  have  the  general  form  given  by  Eqs.  (69) 
and  (70)  .  At  the  receiver  the  received  signal  has  the  form 


H  :  r(t)  =  m  (t)  +  a(t)  tgtO.T] 

o  o 

^ :  r(t)  =  m^(t)  +  a(t)  te[0,TJ 


(73) 


An  exception  to  this  simple  binary  signaling  i:  on-off  keying  and 
the  two  AJ  signaling  schemes  which  consist  of  long  sequences  of  such 
binary  pulses. 
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IV  COHERENT  BINARY  RECEIVER  MODEL 


We  now  consider  how  the  coherent  receiver  processes  the  signal 
r(t),  t e [0,T ]  to  determine  which  binary  signal  was  sent.  If  the  noise 
a(t),  te[0,T]  is  a  Gaussian  process,  then  the  minimum-probabil ity-of- 
error  receiver  has  the  form  shown  in  Fig.  10.  Recall  that  a(t),  te[0,T] 
is  in  our  case  a  non-Gaussian  process  which  is  impulsive  in  character. 


IF  X>0  PICK  H0 
IF  \<0  PICK  H | 

0-7045-4**1 


FIGURE  10  COHERENT  RECEIVER  MODEL— NO  CLIPPING 


Since  the  noise  fluctuation  often  exceeds  the  signal  peak,  in  many  re¬ 
ceivers  a  nonlinear  clipping  function  is  introduced  in  the  receiving 
system  before  the  correlation  stage.  Thus  receivers  typically  have  the 


form  given  by  Fig.  11,  where 

(B 

r'(t)  --Jy(t) 

(“B 


y  ( t )  >  B 

-B  <  y(t)  s  B  .  (74) 

y ( t)  <  -  B 
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On  the  basis  of  our  assumed  noise  model,  the  optimum  maximum  likeli 
hood  receiver  will  certainly  be  different  from  the  above  receiver.  In¬ 
deed,  the  above  receiver  is  not  an  optimum  receiver  for  any  known  noise 
model;  rather,  it  evolved  from  classical  Gaussiar.  noise  receivers  with 
the  nonlinearity  introduced  to  get  rid  of  the  non-Gaussian  impulsiveness 
of  the  noise . 

The  nonlinear  clipping  of  the  received  signal  is  difficult  to 
analyze.  Watt,25  however,  has  shown  that  if  the  clipping  level  is 
greater  than  the  signal  peak  then  most  of  the  slipping  effect  is  on  the 
noise  peaks.  He  has  shown  experimentally  that  the  distribution  of  sig¬ 
nal  and  atmospheric  noise  at  a  level  twice  or  more  above  the  signal  peak 
is  approximately  the  result  of  noise  only.  We  therefore  make  the  assump 
tion  that  when  the  clipping  level  is  above  the  signal  peak  only  the 
noise  signal  peaks  are  clipped.  When  one  considers  the  impulsive  nature 
of  the  noise  and  the  fact  that  the  nonlinear  clipper  is  designed  to  cl  p 
nigh  noise  peaks,  this  assumption  is  reasonable.  The  resulting  signal 
is  shown  in  Fig.  11,  where 


r(t)  =  A  en^  cos^2rrfQt  +  9 ( t ) J  +  m(t) 
r'(t)  =  c ( t )  cosj^-^t  +  9(t)l  +  m(t) 


(75) 


where  m(t)  is  either  mQ(t)  or  m  (t),  and  c(t)  is  given  by  Eq .  (63). 

After  passing  through  the  nonlinear  clipper,  the  received  clipped 
signal  is  now  processed  in  the  usual  correlator  receiver,  which  computes 
a  likelihood  function  on  which  it  bases  a  decision.  Since  correlation 
is  a  linear  operation,  if  the  noise  process 

c(t)  cosTu^f^l  *  0(t)|  (76) 


is  Gaussian,  then  *.  is  a  Gaussiar,  random  variable.  Even  if  the  noise 
process  is  non-Gaussian,  if  the  integration  time  T  is  large  compared  to 
the  correlation  time  of  the  nois,  ,  \  approaches  a  Gaussian 


random  variable  in  distribution.  This  approximation  gets  better  as  T 
increases.  The  correlation  time  of  the  noise  process  is  roughly 


which  for  W  =  1000  Hz  is  0.5  ms.  We  consider  the  signal  time  for  each 
bit  of  approximately  20  ms  to  be  roughly  40  correlation  times.  Recall 
that  for  small  fluctuations  the  atmospheric  noise  process  has  a  Gaussian 
distribution  (Rayleigh  envelope)  and  that  large  fluctuations  which  could 
normally  change  the  distribution  of  the  likelihood  ratio  X  have  been 
clipped  so  that  their  effect  on  X  is  reduced.  These  arguments  lead  us 
to  conclude  that,  although  the  clipped  noise  process  given  by  Expression 
(76)  is  Gaussian  only  for  small  fluctuations,  because  of  clipping  and 
the  large  integration  time,  X  can  be  approximated  as  a  Gaussian  random 
variable.  Certainly  this  approximation  is  less  accurate  for  the  tails 
of  the  distribution  of  X.  Determination  of  the  accuracy  of  this  assump¬ 
tion  is  beyond  the  scope  of  this  study.  We  also  see  that,  because  of 
uncertainties  in  determining  signal  and  noise  parameters,  more  accuracy 
in  the  error  expression  is  not  worth  the  additional  effort. 

From  Fig.  11  we  see  that  X  is  given  by 


T 

j>  r(t)[m° 


(t)  -  mj(t) i  dt 


where 


r*(t)  r  e(t)  eosj^-f^t  +  a(t)j  +  m(t) 


Suppose  that  is  true  and  consequently  that  m(t)  =  m^(t).  Then  we 


T 

/c(t)  cosj^Tri^t  *  (t( t)j  j m^( t )  -  m 


jU)  dt  +  (g  -  1)  PT  .  (80) 


An  error  or  wrong  decision  occurs  f  X  is  greater  than  zero.  This  event 

occurs  with  probabilit\  P  ,  given  by 

c 
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P  =  P  |  \  >  OjK.f 

e  r  <  1> 


Def ining 


we  have 


a'(t)  =  c ( t)  cos^2TTfQt  +  Q(t)j 


T 

a  '  ( t)  jnM  t)  -  m^(t)|  dt 


pe  -  Pr|s  >  (i  -  o)  ptIhJ 


Clearly, 


g  =  0 


2  2 
'g  ~  6 


T  T 

~  "  mi*0,*Hno^  “  n“^2)}  dads  (B6) 


The  noise  process  a'(t)  is  essentially  white  compared  to  the  narrow-band 

signals  m  (t)  and  m,  (t).  Hence 
o  1 

«»> 


is  like  a  Dirac  delta  function  in  .  (86).  This  gives 


— -  T 

2  lUja)  f  r 
K  w  /  [  O 


i  2 

(a)  -  at.  ( .->)  |  da 

1  J 


.2 

K  W 


(1  -  ;)  PT 
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We  know  that 


1  im  2W 
W-® 


sin  2rrWt 
2nWt 


5(t) 


(89) 


For  the  front-ond  bandwidth  W  sufficiently  large  that  nMt)  and  m^t) 
remain  unchanged  by  the  front-end  filter,  this  assumption  is  valid.  But 


Hence 


a,2(r y)  =  c2(t)  cos“^2rrfot  +  £)( t) 


.2r: 


=|c2<t, 


(1  -  ,)  PT 


(90) 


(91) 


The  error  probability  is  then 


(92) 


By  symmetry  we  see  that  the  error  probability  when  m  is  sent  is  the 

o 

same.  Hence  represents  the  total  error  probability.  By  collecting 
all  terms  this  is  given  by 


win?  re 


r>:> 


(93) 


(tM) 


(95) 


2 


y 


and 


10  Jn  log10e 


(96) 


Although  not  all  the  receivers  we  consider  have  a  nonlinear  clipper 
in  the  front,  nonintentional  clipping  often  occurs  when  the  dynamic  range 
of  the  receiver  is  exceeded  by  a  noise  spike.  Consequently,  Eq .  (94) 
should  still  be  valid  without  a  clipper  if  the  clip  level  B  is  chosen  to 
be  as  large  as  the  maximum  signal  level.  Note  that 


and 


1  im  $ 
B-^o 


=  0 


(97) 


hence  for  large  B, 


.  •  n2 

lim  P 
B-a 


^loge(B/A)  \ 

V°'  / 


o 

c“ 


approaches 


=  0 


2  2 
1 im  c  =  A 

B-«> 


3 


(98) 


(99) 


being  the  same  as  the  case  for  no  clipping. 

The  probabil ity-of -error  expression  of  Eq .  (93)  applies  for  coherent 
FSK,  PSK,  and  MSK  signaling  schemes.  We  have  developed  this  expression 
here  mainly  to  present  a  rer>"f,sentat  ive  calculation  of  probability  of 
error  and  to  point  out  the  assumptions  made  throughout  the  remainder  of 
Part  Two.  Although  the  error  calculations  may  differ  for  other  communi¬ 
cation  systems  of  interest,  many  of  the  results  of  this  chapter  are  re¬ 
ferred  to  in  later  chapters.  We  now  consider  specific  modulation- 
demodulation  schemes,  then  compute  error  probabilities,  and  finally  plot 
curves  of  these  probabilities  for  typical  parameters. 
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V  PROBABILITY-OF -ERROR  CALCULATIONS 


This  chapter  develops  probability-of -error  expressions  for  various 
generic  modulation-demodulation  systems.  These  systems  and  the  parameters 
chosen  correspond  roughly  to  existing  and  planned  future  VLF/LF  communi¬ 
cation  systems. 

In  all  cases,  we  consider  modulation-demodulation  systems  for  binary 
signals  denoting  Marks  and  Spaces.  We  denote  T  as  the  time  duration  of 
a  single-channel  binary  signal.  In  some  cases  an  information  bit  may 
require  the  transmission  of  several  channel  bits  (sometimes  referred  to 
as  "chips") .  We  also  assume  that  the  noise  in  the  channel  has  a  con¬ 
stant  spectral  density  ->ver  the  bandwidth,  W,  of  the  receiver's  front- 
end  bandpass  filter.  This  spectral  density  is  denoted  as  Nq  watts  per 
hertz.  In  all  cases  signal  power  is  constant  and  given  by  P  watts  per 
second , 

Throughout  most  of  this  chapter  we  assume  that  the  signaling  time 
is  much  greater  than  1/W,  so  the  test  statistics  upon  which  a  decision 
is  based  can  be  approximated  as  Gaussian  random  variables.  With  clipping 
introduced,  this  assumption  is  even  more  accurate,  since  large  noise 
peaks  are  eliminated,  and  low  noise  levels  which  pass  through  the  clipper 
behave  like  Gaussian  noise.  This  Gaussian  assumption  is  especially  good 
for  an  interesting  case  where  integration  times  are  very  long  compared 
to  1/W .  In  the  Appendix  we  take  the  case  where  T  <  1/W  and  W  is  greater 
than  a  few  hundred  nertz.  HalllF  calls  this  case  the  "short -duration" 
signal  and  notes  that  it  is  very  uncommon  at  VLF  and  LF .  In  this  case 
we  use  the  approach  of  Bello26  to  compute  error  probabilities. 

A .  Binary  On-Off  Keying  (CW) 

On-off  keying  is  perhaps  the  simplest  modulation  technique  for 
transmitting  binary  information.  An  on-off-keyed  radio  system  can  be 
described  as  transmitting  pulses  , 


55 


m(t) 


(100) 


/TP  sin^ti^t  +  $)  for  Mark 

0  for  Space 

where  the  signal  duration  is  T  seconds,  and  the  power  of  the  Mark  pulse 
is  P. 

The  signal  (Mark  or  Space)  is  transmitted  over  a  channel  which  adds 
atmospheric  noise  of  constant  spectral  density  over  the  receiver 
front-end  bandwidth,  W.  We  consider  two  types  of  noncoherent  receiver 
systems,  one  with  clipping  in  the  front  end  and  one  without  clipping. 

We  first  consider  a  noncoherent  receiver  with  no  clipping.  This 
consists  of  a  bandpass  filter  matched  to  the  Mark  pulse  (except  for  an 
unknown  phase),  followed  by  an  envelope  detector  sampled  at  time  T  as 
shown  in  Fig.  12,  where 

h(t)  = 


-\/f 


~  sin  to  (T 
T  o 


t)  te.0,T] 


r(t)  =m(r)  +  o(t) 


h(t) 


1 

ENVELOPE 

DETECTOR 

0-7045-132 


FIGURE  12  NONCOHERENT  DETECTION  OF  CW  PULSE— NO  CLIPPING 


The  sampled  test  statistic,  R  ,  is  given  by 


2  2  2 
R  =  r  +  r 
s  c 


where 


T 

r„  =  -\  /t  /  r(t>  sin  «> 

*  VT  { 


(t)  cos  'jj  t  dt 


(101) 


(102) 
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2 

The  decision  rule  is  to  compare  R  with  a  threshold  as  follows: 

2  2 

R  >  5  choose  Mark 

(103) 

2  2, 

R  i  J  choose  Space 


If  a  Mark  is  sent,  we  have 


and 


r(t)  =  /2F 


sm 


uo  t  +  ft 

l  °  J 


+  a(t) 


(104) 


rg  =  f? T  sin  9  +  ng 


r  =  /PT  cos  9  +  n 
c  c 


(105) 


where 


and 


(106) 


(107) 


are  approximated  as  independent  Gaussian  random  variables  with  zero  mean 

and  variance  N  2 .  This  assumes  that  T  is  large  enough  to  include  many 

independent  samples  of  the  noise.  Under  this  assumption  the  random 
2 

variable  R  is  a  noncentrol  chi-squaro  random  variable  of  order  2.  When 
Mark  is  sent,  the  error  probability  is 
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where 


Q(a,b)  = 


(  ^  1 

t  +  a 

1  exP| - 2 - 


I  (at)  dt 
o 


is  the  Marcum  Q-function ,27  If  a  Space  is  sent,  we  have 


r  =  n 
s  s 


r  =  n 


(109) 


(110) 


and  R  is  a  chi-square  random  variable  of  order  2.  When  Space  is  sent, 
the  error  probability  is  given  by 


I  2  2  2  ) 

P  =P  R=r+r'>  6|Space 

eo  rt  ‘  c  I 


=  exp 


i.iii 

(  Nol 


din 


Assuming  that  Marks  and  Spaces  occur  with  equal  probability,  the  total 
unconditioned  binary  error  probability  is 


I  52) 


exp  -  - 


2  I 


(112) 


Figure  13  gives  this  error  versus  s ignal-to-noise  ratio,  PT/Nq,  for 


various  normalized  thresholds,  b  = 


*  -vK'V 

*  rJ 


Stein  and  Jones25*  show  that  the  threshold  $  that  minimizes  the 
error  probability  is  approximately 


2  v.  FT 
*.  -  N  *  — 

o  4 


(113) 


which  gives  a  minimum  error  probability 


1  1  /  5pf  PT"\  1  |  /  PT  \i 

'  V  *2\>r?expn 


This  is  also  plotted  in  Fig.  13. 


(Ill) 
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SOURCE :  »•'  2«  0-TO49-S7 

FIGURE  13  BINARY  ERROR  PROBABILITY  FOR  ON-OFF 
KEYING-NO  CLIPPING 


For  large  signal -to-noise  ratio 


o 


» 


1 


wo  have 

/  £pt  r  FT  \  ./  fFT  \ 


[3] 


where 


§(a)  = 


/  ^ 


(  ^ ) 

7^  e*p|"  rj dt 


(117) 


is  well-tabulated.21  For  large  signal-to-noise  ratio,  we  thus  have 


P  ^  - 
e  2 


1  J  /PT~  \  1  I  PT  | 

J  (V2No)+2'"P('4Nol 


(118) 


For  the  receiver  with  front-end  clipping,  we  have  the  situation 
shown  in  Fig.  14.  This  receiver  differs  from  the  previous  case  in  that 


\ _ £ 


0-7045-51 


FIGURE  14  NONCOHERENT  DETECTION  OF  CW  PULSE— CLIPPING 


r(t)  is  clipped  for  |r(t)|  >  B.  As  shown  in  the  previous  chapter,  when 

B  2 /TP  ,  (119) 

mcst  of  the  portion  of  the  signal  that  is  clipped  is  the  part  of  the  re¬ 
ceived  signal  that  is  due  to  high  noise  peaks.  We  therefore  assume  that 
the  clipper  output  is  approximated  by 

r'(t)  t-  m(t)  »  a'(t)  ,  (120) 


where 


a 


♦  * 


c(  t ) 


s  in  a.  t 

L  ° 


( 121) 


os  in  Eq .  (79).  The  clipped  noise  is  assumed  to  hove  spectral  density 
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(122) 


N' 

o 


> 


where  c  is  given  by  Eq . 


(94)  . 


Using  the  new  noise  density  after  clipping,  we  have  the  error 
expression 


P 

e 


(123) 


This  expression  depends  on  signal-to-nc ise  ratio,  P/N^,  detection 

threshold,  5^,  clipping  level,  0,  receiver  front  bandwidth,  W,  signal 

duration,  T,  and  noise  parameter,  V..  In  Fig.  15  we  plot  P  versus 

a  e 

P/N  for  optimum  thresholds  and  T  =  0.02,  W  =  1000  Hz,  0=2,  0  =4, 

and  various  values  of  V,  . 

d 

L .  Binary  Frequency-Shift  Keying  (FSK) 

One  of  the  most  common  type  of  binary  signaling  schemes  is  the  FSK 
signaling  scheme  with  a  modulation  index  of  one.  This  scheme  results 
in  a  Mark -Space  signal  sequence  which  has  constant  amplitude  and  con¬ 
tinuous  phase.  The  two  signal  pulses  are 


( t) 

=  /’ST5  sin 

L  V  *  vxt] 

for 

Mark 

m  (t) 

=  /2T*  sinl 

it-  ixt 

for 

Space 

o 

o  J 

(124) 


where 

a.  is  the  center  Ireuuuncv  and 
o 

vt  =  "  T.  both  in  radians  per  second. 


This  frequency  deviation,  _)i,  corresponds  to  a  modulation  index  of 

one  and  results  in  the  two  pulses  m,(t)  and  a  (t)  being  orthogonal. 

I  o 

That  is. 


I  m  ( t  i  is ,  ( t )  dt  .  0 
I  o  I 

■4) 


(125) 


j 

i 

i . 


t 


| 


i 


I  * 
i 

*  ■ 
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BINARY  ERROR  PROBABILITY,  p,  binary  error  probability, 


FIGURE  15  BINARY  ERROR  PROBABILITY  FOR  ON-OFF 
KEYING — CLIPPING 
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ye  consider  coherent  and  noncoherent  detection  of  these  FSK  signals 
under  the  two  conditions  of  clipping  and  no  clipping  xn  the  front  end  of 
ihe  receiving  system.  Classical  coherent  detection  without  any  nonlinear 
clipping  is  diagrammed  in  Fig.  16,  where 


(126) 


0-7045-51 

FIGURE  16  COHERENT  DETECTION  OF  FSK— NO  CUPPING 


The  sampled  tost  statistic  is  given  by 


ri  ‘  ro  ' 


where 


T 

ri  /  *■(*)  yf -su{v  *  ut]  d* 

ro  f  r(°  yf  sin[Aoi  *  -  it] dt 

Jo 


(  127, 


The  decision  rule  is  to  compare  X  with  zero,  as  follows, 


X  >  0  choose  Mark 


X  0  choose  Space 


(128) 


where,  as  usual,  we  assume  that  Marks  and  Spaces  occur  with  equal  proba¬ 
bility.  Note  that  the  coherent  receiver  is  merely  another  representa¬ 
tion  of  the  correlation  receiver  of  Fig.  10.  The  error  probability  with 

2  2  2  an 

no  clipping  is  given  by  Eq .  (93)  with  p  =  0  and  c  =  A  e  .  Hence 


P 

e 


(129) 


since 


9  20n 

w  n 

A  e 


=  2N  W 
o 


(130) 


With  clipping  in  the  front  end  we  I.ave,  from  Eq .  (93)  with  p  =  0, 


P 

e 


(131) 


where 


Vd  =  10  Jn  log10C 


(132) 


(133) 


Choosing 


B  = 


(  13-1) 
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we 


rewrite  in  terms  of  parameters  P y  g,  T,  and  W: 

i ^ 

PT 


Pe  =  f| 


"ol1  -  i<D)  +  i>2  iTw  5<E>| 

1  O 


(135) 


where 


20  log  (0//W)  +  10  log  (P/N  )  -  2V  log  10 


(136) 


20  n  /  -  loge  10 


and 


20  log  ((3//W)  +  10  log  J  P/N  )  +  2V  log  10 
E  =  - 2 - ® -  .  (137) 

20  iose 10 

Figure  17  plots  P^  for  FSK  signals  for  various  parameters  of  interest. 
Here  3  represents  the  clipping  lever  a  ove  the  signal  peak  of  /TP.  For 
example,  0=2  met  ns  tha*  the  clippinr  level  is  set  at  B  =  2/2P  or  twice 
the  signal  peak  amplitude. 

For  p'  ncoherent  detection  without  clipping,  we  have  filters  matched 
to  the  Mr*  a:d  pace  pulses  (except  for  an  unknown  phase),  followed  by 
envelope  detectors  sampled  at  time  T,  as  shown  in  Fig.  18,  where  the 
sampled  test  statistic  is  given  by 


•> 


2 


Is  lc 


2  2  2 

R  =  r  *  r 
o  os  oc 


(138) 


The  terms  r  and  r  are  the  sine  and  cosine  components  of  the  received 
Isle 

signal  relative  to  the  Mark  frequency.  Similarly,  r  and  r  are  the 

os  oe 

components  relative  to  the  Space  frequency. 

Recalling  that  m,(t)  and  m  (t)  are  orthogonal  or  uncorrel atod ,  we 
1  o 

see  that  the  matched  filters  h  (t)  and  h  (t)  arc  also  orthogonal.  This 

1  o 
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h,(t) 


ENVELOPE 

DETECTOR 


hQ(t ) 


ENVELOPE 

DETECTOR 


I - -  I  I  I  0-7048-47 

FIGURE  18  NONCOHERENT  DETECTION  OF  FSK— NO  CLIPPING 


means  that  the  outputs  at  time  T  of  the  matched  filters  have  noise  com¬ 
ponents  that  are  uncorrelated  and,  since  they  are  assumed  to  be  Gaussian 

random  variables,  independent.  Assuming  a  Mark  is  sent,  we  then  have 

2  2 
R  ,  a  noncentral  chi-sauare  random  variable,  and  R  ,  a  central  chi- 
1'  '  ’  o* 

square  random  variable.  Perhaps  a  conceptually  simpler  representation 
of  the  noncoherent  receiver  process  is  that  given  in  Fig.  19,  where 


-1  f 

r,c 

Vv 

Jo 

♦  R 


A  . 

r 

r.s 

jQv 

J* 

♦) — ■‘A 


•T  r  ^  r  1 


4>  in 

o* 


*—  0-7048-48 

FIGURE  19  NONCOHERENT  FSK  DETECTOR  REPRESENTATION 


(139) 


1c 


Is 
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OS 


(t)  =  -v 

/§  cosU  t 

+  Atut 

V 

'  T  [  o 

(t)  = 

/ I  sin joj  t 

+  Aiut 

A/ 

'  T  [_  o 

(t)  = 

-  COS  ID  t 

-  4uut 

V 

!  t  L  ° 

(t)  = 

/-  sin U)  t 

-  Au)t 

'  T  [  o 

If  a  Mark  is  sent,  the  probability  of  error  is  found  from  the  dis- 
2  2 

tribution  of  R,  -  R  .  First,  define  vectors 
1  o  ’ 


(140) 


Since  a  Mark  is  sent,  consists  of  signal  plus  noise,  while  r^  consists 
of  noise  alone.  Hence, 


r. 

/PT  cos  q  +  n. 

lc 

lc 

Jls_ 

-  /PT  sin  9  +  n^s 

m  +  n^ 


r 

n 

oc 

oc 

r 

n 

os 

os 

(142) 


where  ft  is  an  unknown  phase  term,  and  n,  ,  n,  ,  n  ,  n  are  independent, 

°  ’  lc'  Is’  oc’  os  ' 

zero-mean  Gaussian  random  variables  with  common  variance  N  >2.  Using  the 

o 

notation 


,,,22  2 
r  =  r,  ♦  «\ 

1  — 1 1 1  lc  Is 


y 
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we  have 


*1  =  Hiill2  -  He  +  iill2 


<-  HiJI2  -  llJ„ll2 


The  probability  of  error  is 


p  =  P  A  <  0  H 
e  r|  1| 


P  -  R2  <  0|H  j 
r  1  o  1  1) 


pr 1 1 1e  +  Ex 1 12  -  llEoll2 «  °| 


pr!ll“  -  “Jl2  <  llsoll2! 


[  1  pr!iiSoU2  >  Hi  *  — 1 1 ' z i  — i 


LL 


|  — 1  =  d2 


But 


t  .  ,  .  .2  2t  |  R  | 

pr  1 1  iHol  1  *  R  |  =  exp|-  iT|  ’ 


SO 


P  r  J I S  Do  H  2  ?  i!m  +  n«l|2In.  =  or 


i-  -i"  i-i  - 


»"  I 


„  i_  l la  +  al 

xPi  -LJ — jj — 


Also, 


Pn  <2>  a  —  exP 


2-1 


1  \  IWII2^ 

—  exp  -  IL-ii- 
o  l  of 
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(143) 


(144) 


(145) 


(146) 


(147) 


Hence 


or 


(149) 


Although  this  expression  is  derived  for  FSK  signals,  it  applies  in 
general  for  noncoherent  detection  of  binary  orthogonal  signals, 

With  nonlinear  clipping,  we  can  follow  arguments  similar  to  the 
earlier  results  ^o  get  an  error  probability 

1  I  PT  ) 

Pe  =2  eXP(-2SJ)  '  (150) 

where  N'  satisfies 
o 

N'W  =  ^  c2  ,  (151) 

o  A 


or,  in  terms  of  basic  parameters, 

K  -  NoJ‘  -  ,<D>  *  92(iT«),lE)j  '  (152> 

where  D  and  E  are  given  Eqs.  (136)  and  (137).  Equation  (150)  is 
plotted  for  various  parameters  in  Fig.  20. 

C .  Binary  Phase-Shift  Keying  (PSK) 

Although  PSK  is  not  commonly  used  at  VLF/LF,  we  present  it  here 
because  it  is  similar  to  the  minimum-shift -keying  (MSK)  scheme  now 
planned  for  several  new  VLF/LF  systems.  For  binary  PSK,  the  two  signals 
are 
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(153) 


m^d)  =  /27  sin  uo^t  for 


Mark 


m  (t)  =  -  /2F  sin  tu  t  for  Space 
o  o 


Note  that  for  this  signal  set  we  have  simply 


P  =  -  1 


(154) 


Coherent  detection  without  any  nonlinear  clipping  is  simply  represented 

in  the  receiver  of  Fig.  10,  and  the  error  probability  is  given  by  Eq  . 

2  2  2  2 

(93)  with  p  =  -  1  and  c  =  A  e  an .  Hence, 


(155) 


Since 


2  < 

A  e  =  2N  W 
o 


(156) 


We  plot  this  curve  in  Fig.  21  where 


(157) 


With  clipping  in  the  front  end,  we  have  from  Eq.  (93)  with  0  =  -  1, 


P  =  4 
e 


(V-?)  ■ 


(158) 


where,  as  before. 


2  2 
e  =  A  c  1  -  i 


2  nl,  .1"**'*'"  -2<\ | 
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(159) 
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FIGURE  21  BINARY  ERROR  PROBABILITY  FOR  COHERENT  FSK — 
NO  CLIPPING 


Again,  choosing  the  clipping  level, 

B  =  «  /TTP’  (160) 


we  can  write  P  in  terms  of 


parameters 


P 

o 


rfl 


(161) 


where  l)  and  E  are  defined  in  Eqs .  (136)  and  (137). 

Noncoherent  detection  of  PSK  signals  does  not  exist,  except  perhaps 
-  ith  a  phase  reference  signal  .  Suppose,  for  example,  beside  the  Mark- 
Space  signals  we  liave  a  reference  signal  given  by 


q { t )  ' JtP  sin  j.  t 

o 


(162) 


| 
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■  ■  IJWe* 


which  is  first  transmitted  during  time 

[-T,  0  ]  .  ' 

Then,  during 

time  interval 

[-T,0], 

r'(t)  =  /12P  sin(  ujQt  +  9  j 

+  n'(t) 

t«[-T,0] 

(163) 

is  received,  while  during  [0,T], 

r(t)  =  /PP  sin|uoot  +  9)  + 

n(  t) 

for  Mark 

(164) 

r(t)  =  -  /2F  sin(iuot  +  9) 

+  n(t) 

for  Space 

is  received. 

Thi  situation  is  equivalent  to  receiving  a 

total  signal 

(r'(t) 

t  e [-T,o  ] 

rTU)  =| 

(165) 

/  r(t) 

t  €  [0 ,  T 1 

in  f-T,T]  on  which  to  base  a  decision 

on  whether 

a  Mark  or  a 

Space  has 

been  sent.  Defining  signals 

(qU) 

te[-T,0] 

n1  < t)  *] 

(166) 

t  6  [0  .  T  ] 

and 

iQ(t) 

tC[‘T,0  j 

mT(t)  =| 

9 

(167) 

0  ) 

I  m  (t) 

t  c  [0  ,  T  1 
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wo  have 


my(t)m^(t)dt 
1  o 


I  q2(t)dt  + 

J-T 


m  (t)m, ( t)dt 
o  1 


(168) 


=  PT  -  PT 


=  0 


We  have  the  situation  where  two  orthogonal  signals  are  transmitted  over 
[-T,Tj  with  a  random  phase  9.  Noncoherent  detection  of  such  signals  is 
precisely  the  same  as  noncoherent  detection  of  FSK  signals  (orthogonal 
signals)  of  total  duration  2T.  Hence,  from  Eq  .  (149),  we  have 

1  (  PT  I 

P*  '  2  “P|  iT|  U69) 


for  noncoherent  detection  of  PSK  signals  with  a  reference  signal  as 
stated  above.  For  the  case  with  clipping  we  have,  from  the  results  of 
FSK,  the  error  probability 


P 

e 


y 


where 


S' 

o 


°l 


J(L>) 


(170) 


(171) 


The  most  common  reference  signal  is  the  previously  transmitted  Mark  or 
Space  signal.  This  means  that  each  transmitted  signal  is  modulated  rela¬ 
tive  to  the  previously  transmitted  signal.  This  efficient  mode  of  non¬ 
coherent  operation,  called  d i f ferent iol -phase  shift  keying  (DPSK) , 
results  in  the  binary  error  probabilities  given  above. 
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D  .  Binary  Minimum-Shift  Keying  (MSK) 

Minimum  shift  keying  is  a  patented  scheme  developed  by  the  Collins 
Radio  Company.  This  modulation  scheme  gives  a  signal  in  the  channel 
that  looks  like  FSK  with  a  modulation  ;ndex  of  m  -  1/2.  Normally  FSK 
modulation  with  m  =  1/2  would  require  discontinuous  phase  changes,  which 
would  be  difficult  for  VLF/LF  transmitters  to  handle.  The  MSK  scheme 
gives  continuous  phase  and  an  improved  performance  equivalent  to  that  of 
PSK  modulation. 

Consider  a  basic  information  source  that  generates  a  Mark  or  Space 
once  every  T  seconds.  This  binary  information  sequence  is  alternatively 
used  to  modulate  two  PSK  modulators  that  art-  *r/2  radians  out  of  phase 
with  each  other.  This  results  in  two  PSK  signal  sequences  where  phase 
changes  occur  at  odd  times  (T,  3T,  5T,  . . .)  in  one  and  at  even  times 

(0,  2T,  4T,  ...)  in  the  other.  The  tv e  PSK  signals  are  amplitude  modu¬ 

lated  by  /2P  sin  uo^t  and  /2T>  cos  x^t,  respectively,  and  summed  before 
transmission.  See  Fig.  22,  where  =  n/2T, 


FIGURE  22  MSK  MODULATOR  REPRESENTATION 
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,  s  in 


jl,  t  sin  1 1  t 
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a  I  +  cos  x,  t  cos:  x  t 
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(172) 
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The  signal  S(t)  shown  above  appears  as  a  continuous  phase  FSK 
signal  with  a  modulation  index  of  m  =  1/2.  Coherent  detection  is  accom¬ 
plished  by  separating  the  signal  plu3  noise  into  two  separate  PSK-like 
signal-plus-noise  sequences  which  are  detected  separately  as  shown  in 
Fig.  23,  where 


*s(t) 

<tc(t) 


"\/f  si“  “i 


t  sin  u)  t 
o 


t  COS  u)  t 

o 


(173) 


and 


®s(t)  <tc ( t)  dt 


^  0 


it  ( t)  e  ( t)  dt  =  0 
s  c 


(174) 


FIGURE  23  COHERENT  MSK  DETECTOR— NO  CLIPPING 


Since  each  side  of  the  detector  is  essentially  a  PSK  detector,  the 
performance  is  that  of  coherent  PSK  detection  with  average  power  1/2  P 
and  signal  duration  2T  which  results  in  total  signal  energy  PT  per  trans¬ 
mitted  Mark  or  Space.  Hence,  for  each  transmitted  bit,  the  error  proba¬ 
bility  for  this  coherent  detector  is 
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when  no  clipping  is  employed.  This  is  plotted  in  Fig.  24  for  T  = 


FIGURE  24  BINARY  ERROR  PROBABILITY  FOR  COHERENT  MSK — 
NO  CLIPPING 


With  clipping  we  have 


P 

e 


2PT _ 

*(D)  + 


} 


just  as  in  normal  PSK  signaling  schemes.  These  probabilities  ore 
in  Fig.  25. 


(175) 

0.02. 


(176) 

plot  ted 


We  were  able  to  treat  each  side  of  the  detector  separately  because 
of  the  orthogonality  property  of  «s(t)  and  ©c(t)  shown  in  Eq .  (171)  and 
the  fact  that  wo  have  exact  knowledge  of  the  signal  phase.  We  consider 


now  the  difficult  problem  of  operating  in  a  differential  mode  similar  to 
differential  phase  shift  keying. 

Suppose  we  have  a  reference  signal 


q(t)  =  rZ$  J  sin  (iLt  sin  uj  t  +  cos  ai,t  cosltu  t  +  0,)J  (177) 

)  1  o  1  'o  y )  | 


which  is  transmitted  during  [-T,Tl  and  received  as 


r'(t)  =  /2P  sin  uo^t  sinjx^t  -  g)  +  cos  u^t  cos(u>ot  +  Q1  -  9)  +  n'(t) 


Def ining 


(178) 


«  (t) 

ss 


4>  (t) 

sc 


<5  (t) 

CS 


0  (t) 

cc 


t  sin  u)  t 
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•  y^in  ”1 

-  yi" sin 

=  "V? cos  “1 

[2 

=  ~\ij  c°s  c°s  uj^t 


t  cos  UJ  t 
o 


t  sin  x  t 
o 


(179) 


we  now  consider  how  uncertainty  in  the  phase  of  the  transmitted  signal 
can  cause  some  mutual  interference  between  the  two  PSK  sides  of  the 
detector.  Consider 
T 

(180) 


•  L " 


(t)«  (t)  dt 

ss 


T 

j  sin2  *ix  sin[xGt  "  s] 


sin  x  t  dt 
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/-T 

T 


+  sin  tjt  cos  x^t  cos[xot  *  9'  -  qj  sin  jM  di 
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Similarly, 


r'  =  -  /PT  sin  9  +  n' 
c  c 


+ 


0 


2 

-  /FT  cos  9 


-  /FT  cos  9 

TT  ' 


0'  =  0  in  [-2T.2T] 

y 

0'  =  tt  in  [-2T , 2T ] 

y 

S'  =  0  in  [-2T,  0],  0' 
0y  =  tt  in  [-2T,  0] ,  0^ 


tt  in  [0, 2T] 

0  in  [0, 2T] 


(184) 


Here  n'  and  n'  are  independent  zero-mean  Gaussian  random  variables  with 

variance  N  /2 .  Defining  the  random  variable  that  takes  on  one  of  three 
o 

values  with  the  following  probabilities, 


6 


(185) 


we  can  represent  the  sine  and  cosine  components  of  our  reference  signal 
as 


_  2  A 

r'  =  /FT  cos  0  +  n#  +  — TFT  sin  0 

S'  S  TT 

_  2  f) 

r'  =  -  /FT  sin  Q  +  n'  +  —  /TT  cos  9 


(186) 


This  shows  that,  besides  the  noise  terms  n'  and  n^,  the  reference  signal 
can  be  corrupt*. '  by  a  strong  interference  term  due  to  the  signals  in  the 
other  half  of  the  detector.  This  means  that  using  MSK  on  a  bit-by  bit 
differential  mode  is  impractical.  We  shall  see,  however,  that  taking  a 
long  stream  of  such  channel  bits  (chips)  to  form  a  reference  for  another 
long  stream  of  bits  (or  chips)  can  result  in  good  performance. 
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Suppose  that  we  now  take  N  binary  MSK  channel  bits  or  chips  to 
transmit  a  single  Mark  or  Space.  This  can  be  accomplished,  for  example 
by  transmitting  a  pseudo-random  binary  sequence  of  chips  for  Mark  and 
its  negative  for  Space.  Assume  that  we  take  N  pseudo-randomly  generated 
binary  MSK  chips  to  transmit  a  Mark  or  Space,  and  suppose  that  the  N 
bits  of  the  previously  transmitted  Mark  or  Space  are  used  as  a  phase 
reference.  This  means  that  we  have  a  noncoherent  differential  mode  of 
operation  on  N  chips,  similar  to  the  differential  PSK  scheme  described 
earlier . 


Collecting  the  N  sine  and  cosine  components  of  our  previously  trans¬ 
mitted  signal,  we  have 


N 


k=l 


=  N/FT  cos  q 


N 

E  ns<k> 

k=l 


(187; 


and 


♦s'" sl"  9E  5k 

k=l 


R'  -  V'  r' 
c  c 
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k=l 
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9  *E  nc' 


N/FT  sin  0  +  )  n'(k) 

k=l 
N 


+  -  'FT  sin  n  )  6, 

TT  Z _ ,  k 


E 

k=l 


(188) 


as  the  sine  and  cosine  components  of  our  reference  sequence  where  the 
index  k  corresponds  to  the  kth  chip  in  the  reference  sequence.  The  terms 
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and 


N 


k=l 


(189) 


(190) 


are  Independent  zero-mean  Gaussian  random  variables  with  variance  NNq/2 . 
Assuming  that  chip  sequences  are  pseudo-randomly^  generated,  we  can  take 
the  random  variables  {6^}^  ^  to  be  independent,  so  that  the  sequence 


D' 


(191) 


approaches  a  Gaussian  random  variable  in  distribution  as  N  gets  large. 
The  sine  and  cosine  components  of  the  reference  signal  are  thus 


R'  =  tyPT  cos  0  +  N'  +  -  /FF  sin  0  D; 

S  S  TT 

2 

r'  =  -  N/PT  sin  Q  +  N'  +  -  PT  cos  9  D' 

c  C  TT  ■ 


(192) 


where  N',  N^,  and  D'  are  independent  zero-mean  Gaussian  random  variables 
(D#  being  assumed  to  be  Gaussian)  with  variances 


N'  2  =  N'  2  =  NN  /2 

SCO 
D2  =  N/2 


(193) 


Note  that  the  noise  components  due  to  the  interference, 

9 

-  -TT  sin  a  D' 

TT 

2  , 

-  /pt  cos  a  D 


( 19-*) 


8-1 


form  a  vector  that,  is  orthogonal  to  the  signal  vector  formed  by 


t 


a.*  io‘2  L- 


2  »o* 

o 

cc 

a 


oc 

< 

£  i0'* 


-20  -15  -1C  -5  0  5  10  15  20 

SICNAU-TO-NOISE  RATIO.  P/N,  - dB 

0  0-T045 

FIGURE  26  BINARY  ERROR  PROBABILITY  FOR  DIFFERENTIAL  MSK 


where,  as  before, 


N'  =  N  a  -  $(D)  +  0 
o  o  I 


(202) 


M-ary  Minimum  Shift  Keyit 


We  now  consider  schemes  that  take  N  binary  MSK  chips  to  transmit 
one  of  M  messages.  In  general,  the  performance  of  these  schemes  is 
difficult  to  evaluate,  so  we  must  resort  to  finding  upper  bounds  on  the 
error  probabilities.  In  this  section,  we  restrict  ourselves  to  two 
M-ary  signal  schemes,  coherent  orthogonal  signals,  and  coherent  pseudo- 
randomly  generated  signals. 

Suppose  that  the  M  signals  consisting  of  a  sequence  of  binary  MSK 
chips  are  represented  by 
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N-l 


s±( t)  =y~|  stJp(t  -  jT)  t«[0,NTl  i  =  1,  ...,  M  ,  (203) 


where 


j=0 


I  1  ^ 

5u  *i-i’  and 


p(t)  Is  a  pulse  of  duration  T  and  energy  PT. 


Also,  let 


r(t)  te[0,NT] 


be  the  received  signal  which  consists  of  the  transmitted  signal  plus 
add  it iv 
f ining 


additive  atmospheric  noise  with  spectral  density  Nq/2  as  before.  De- 


NT 

I'll  2-f  * 

Jq 


(t)  dt 


(204) 


for  any  time  function  f ( 1 )  over  [0,NT],  we  consider  events 


e. 


ik 


|l!r  ~  skl!2  <  Hr  -  si 1 12  si  is  sentj 


(205) 


,  th 


for  i,  k  =  1,  2,  M.  If  the  i  signal  is  transmitted  and  an  optimum 

white  Gaussian  noise  receiver  is  used  (as  in  the  case  of  interest),  the 

kth  signal  is  mistakenly  assumed  to  have  been  sent  by  the  receiver  when 

the  event  f  occurs.  This  event  occurs  with  probability 
ik 

,.e 

rrik*  “  rj "  >  2  \ 


■ 


(206) 


th 


In  general,  when  the  i  signal  is  transmitted,  an  error  occurs  if  one 
or  more  of  the  events 


t. 


ik 


k  ,  1,  2.  . . M 


k  *  i 


(207) 
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I 


occurs.  Hence,  the  error  probability  when  the  i 
is  given  by 

M 

,  .  )  V 

P  £ls.  is  sent 


th 


rr <  i 


'! 


■ p*5J s 


k=l 

Mi 


A  well-known  property  of  probabilities  is  that  the 
of  events  is  less  than  the  sum  of  the  probabilities 
we  have  a  bound 
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Mi 


r(  ikl 
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When  there  is  an  orthogonal  signal  set,  then 
NT 


l 


s  .  ( t)  s,  (t>  dt  =  0  k  t  i 
l  k 


and 


|si  "  sk||  =  2NPT  k  f  i 


Thus,  for  the  orthogonal  signal  set,  we  have 
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P  £|s,  is  sent  5 
n  1  i  I 


5<vn 


Mi 


or 


P^. }  £  1 s  ^  is  sentj  < 


■•(V?) 


signal  is  transmitted 

(208) 

probability  of  a  union 
of  each  event.  Hence, 

(209) 


(210) 


(211) 


(212) 


(213) 


(214) 
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Since  this  bound  does  not  depend  on  which  signal  is  sent,  the  error 
probability  when  any  signal  is  sent  has  the  same  bound.  The  probability 


oi  error,  P  ,  when  any  signal  is  sent  is  thus  bounded  by, 


which  is  plotted  in  Fig.  27. 


(215) 


FIGURE  27  CHARACTER  ERROR  PROBABILITY 
OF  MSK  SEQUENCES 


When  the  signals  arc  not  orthogonal  but  are  generated  pseudo- 
rnndomly  by  independent  pseudo-random-sequence  generators,  Eq .  (213) 
still  applies.  Hence,  when  the  i*h  signal  is  sent,  the  error  probability 
is  bounded  by 


P  £  s.  is  sent  £ 
r  i 


gf^) 


(216) 


Mi 


In  general,  we  are  interested  not  in  the  performance  of  a  given  pseudo- 
randomly  generated  signal  set  but  in  the  average  performance  when  a  large 
class  of  such  signal  sets  are  used.  Hence,  when  s^  is  sent,  we  are 
interested  in  the  average  error  probability  where  we  average  over  the 
whole  ensemble  of  pseudo-randomly  generated  signal  sets.  Taking  averages 
of  Eq .  (216)  we  have, 

Pr|e|s.  77  sent}  ♦(— Tfr--)  ■ 

k=l  '  °  ' 

Mi 


(217) 


The  error  function  has  a  bound 


$(x)  <  exp 


I.  <2I 


|  2  |  ’ 

so  that  Eq.  (217)  can  be  further  bounded  by 
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Note  that 
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We  assume  that  the  ensemble  of  signal  sets  is  generated  in  such  a  manner 
that  the  components  of  any  signal  are  essentially  independent  of  each 
other  and  of  other  signal  components  and  take  on  either  +1  or  -1  with 
equal  probability.  Hence, 


2  + 


exp 


L  nl 

I  »ol 


(222) 


The  total  probability  of  error,  given  that  s^  is  transmitted,  is  bounded 

by 


£  s. 
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or 
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1 

2  GX 
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(224) 


is  the  bound  on  the  probability  that  any  signal  transmitted  gives  an 
error.  Equation  (224)  can  be  written  in  a  more  familiar  form  as, 

P  M  exp  N  log 

•  I 

This  is  plotted  in  Fig,  27. 

In  both  Eq .  (215)  and  Eq .  (225),  for  the  case  of  front  end  non¬ 
linear  clipping,  the  term  is  again  replaced  by 

k  ■  "o)1  -  *iD>  *  32(sfs-)i<E>j  i226> 


as  before. 
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In  the  preceding  derivation  of  the  performance  of  typical  equip¬ 
ment  modems  in  the  presence  of  atmospheric  noise,  improvement  is  shown 
to  result  from  wideband  clipping  at  the  front  end  of  the  receiver.  The 
development  accounts  for  the  reduction  in  impulsive  noise  by  the  clipper 
and  hence  predicts  an  improvement  in  signal-to-noise  ratio  in  the  stages 
of  the  receiver  following  the  clipper.  While  the  predicted  improvement 
is  valid  when  the  clipping  level  is  well  above  the  rms  noise  level  in 
the  wideband  portion  of  the  receiver,  when  the  noise  clipper  is  active 
an  appreciable  fraction  of  the  time  these  predictions  become  inaccurate. 
One  reason  for  this  inaccuracy  is  that  account  is  not  taken  of  the  loss 
in  signal  energy  as  a  result  of  the  operation  of  the  noise  clipper. 

An  approximate  correction  may  be  made  as  follows: 

let  P  =  probability  that  the  noise  level  exceeds  the  clipping 
B 

level,  B. 

Then,  during  the  time  that  the  clipper  is  saturated  by  large  noise  pulses, 
no  signal  power  is  transmitted  through  the  clipper,  and  the  loss  in  sig¬ 
nal  power  is  given  by: 

P_  1 
P'  =  1  -  PB 

All  appropriate  curves  in  this  report  include  this  correction. 
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VI  CONCLUSION 


We  have  computed  the  performance  of  several  generic  VLF/LF 
modulation-demodulation  schemes,  using  the  atmospheric  noise  model  de¬ 
veloped  in  Part  One.  In  general,  these  calculations  assume  th£t  the 
integration  time  T  (baud  time)  is  large  enough  so  that  the  integrator 
output  is  Gaussian  in  distribution.  When  clipping  is  performed  before 
integration,  this  assumption  is  even  better.  For  the  case  of  no  clipping, 
this  analysis  results  in  classical  white  Gaussian  noise  error  proba¬ 
bilities.  With  clipping,  the  detail  properties  of  atmospheric  noise 
! 

become  important,  and  such  parameters  as  front-end  bandwidth  (W) ,  clip 
level  relative  to  signal  peak  (g),  and  noise  voltage  deviation  (V^)  be¬ 
come  important . 


93 


APPENDIX 


Here  we  consider  coherent  PSK  and  MSK  performance  when  the  inte¬ 
gration  time  T  is  small  compared  to  the  atmospheric  noise  correlation 
time  and  when  no  clipping  is  introduced.  First  consider  the  PSK  signals 


S  (t)  ^  /"ZF  sin  a]  t  0  t  5  T 

o  1  o 


S,(t)  =  -  sin  'jo  t  0  £  t  £  T 
1  o 


(227) 


Suppose  that  the  signal  S  (t)  is  transmitted.  The  received  signal  is 

o 

then 


.  n(t) 

r(t)  =  A  e  cos 


[v  +  e]  +  So(t) 


(228) 


The  coherent  receiver  correlates  the  received  signal  with 


AS(t)  =  S  (t)  -  S,(t) 
o  1 
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to  obtain  the  test  statistics 
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We  now  make  the  assumption  that  T  is  small  enough  so  that  n(t)  and 
6  are  essentially  constant  over  the  time  interval  [0,T].  Therefore 


K  =  2PT  +  2A  e  JTP 


f 

I  cos  uu  t  +  9  sin  w  1 

/  L  0  J  ° 

Jo 


(232) 


=  2PT  +  JTP  A  e  T  sin  0 


An  error  occurs  if  \  £  0  when  S  (t)  is  transmitted.  Hence  the 

o 

error  probability  is 


P  =  P  \  £  0  S  is  sent 
e  r  o  ! 


=  P  !  2PT  +  JTP  Ae  T  sin  9  £  oj 


(233) 


„  |  n  .  JTP\ 

=  Prje  sm  9  £  -  i—  j 

Since  we  assume  that  9  is  completely  unknown,  we  assume  that  it  has 
uniform  probability  over  [0,2tt].  Consequently,  its  probability  density 


p9<«> 


—  0  £  a  £  2tt 


elsewhere 


(234) 


Since  sin  9  can  be  positive  or  negative  with  equal  probability,  by 


symmetry  we  have 


P  =  P  !en  sin  9  2 
e  rj  A  J 


(235) 


Clearly,  when  n  £  9  £  2n, 


e  sin  9  >  *— — 


(236) 


96 


(237) 


TT 


Ry  symmetry  v;e  note  that  the  same  error  probability  applies  when 
S^(t)  is  transmitted.  Hence  the  total  binary  error  probability  for  co¬ 
herent  binary  PSK  with  integration  time  T  smaller  than  the  noise  corre¬ 
lation  time  is  given  by 


P  =  — 


1  f  .heir^c, 

E 1  \i^ 


d0 


(238) 


This  error  probability  is  generally  difficult  to  evaluate  except  by 
approximation  and  computer  techniques.  Coherent  MSK  also  as  this  form 
for  the  probability  of  making  a  binary  error. 
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havior  of  atmospheric  noise  becomes  more  important,  and  results  depend  cn  the 
detailed  structure  of  the  noise. 
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